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Abstract 



PLh ' Smoothed Wigner transforms have been used in signal processing, as a regularized 

■^ . version of the Wigner transform, and have been proposed as an alternative to it in the 

(-^ I homogenization and / or semiclassical limits of wave equations. 

We derive explicit, closed formulations for the coarse-scale representation of the action 
of pseudo differential operators. The resulting "smoothed operators" are in general of 
infinite order. The formulation of an appropriate framework, resembling the Gelfand- 
Shilov spaces, is necessary. 
^ ' Similarly we treat the "smoothed Wigner calculus". In particular this allows us to 

0^ , reformulate any linear equation, as well as certain nonlinear ones (e.g. Hartree and cubic 

pj I non- linear Schrodinger), as coarse-scale phase-space equations (e.g. smoothed Vlasov), 

f^ ■ with spatial and spectral resolutions controlled by two free parameters. Finally, it is 

T-h . seen that the smoothed Wigner calculus can be approximated, uniformly on phase-space, 

CD ' by differential operators in the semiclassical regime. This improves the respective weak- 

^i: ■ topology approximation result for the Wigner calculus. 
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1 Introduction and statement of the main results 

Homogenization is an increasingly important and diverse paradigm of applied mathematics. 
It is fair to say that it often consists in the reduction of a "complicated" problem, typically 
involving multiple scales, to an effective problem which describes correctly certain coarse-scale 
features, while appropriately averaging "less important" ones, without having to keep track of 
them explicitly. The simpler, effective problem is then more amenable to analytical and / or 
numerical treatment. 

The Wigner transform (WT) has been used extensively in the homogenization of wave 
problems, and notably in semiclassical asymptotics for the Schrodinger equation. As we briefly 
mention in the abstract, this is where the motivation for the smoothed Wigner transform 
(SWT) and for the study of "smoothed operators" comes from. However, the development of 
the "smoothed calculus" is sophisticated enough on its own, and completely independent from 
the specifics of the WT or the SWT on the technical level. 

Because of that, and in order to make the presentation friendlier to the reader, the paper 
and its introduction have two parts. First we discuss the smoothed calculus: in Section [LD we 
motivate the result, and outline its simplest formulation. In Section [L2] we briefly outline the 
application of the smoothed calculus to a SWT-based scheme for the reformulation of general 
classes of wave equations to coarse-scale kinetic problems in phase-space, i.e. in a space of 
position X and momentum / wavenumber k. 



1.1 Deriving smoothed equations 

On the technical level, the derivation of equations for the SWT consists in "smoothing" the 
well known Wigner equations. That is, starting from a well-defined system of equation^ 

Lw = (1) 

for an appropriate function w, we want to derive the equations governing a smoothed version 
of w, symbolically ^w. In other words - and more generally - we want to smooth equation ([T]) 
and commute correctly the smoothing with the operator, 

Lw = ^ ^Lw = ^ L{^w) = 0, (2) 

so as to get a closed problem for the smoothed function w = ^w. As we will see in the sequel, 
building a useful, practical theory using the "smoothed" problem 

Lw = 0, (3) 

involves more work than just deriving it. 

In the case of smoothed Wigner transforms, which is our concrete motivation, it must 
be emphasized that there are sound mathematical and physical reasons to believe that the 
smoothed equations are useful (at least for certain problems). In addition, computational 
aspects (in particular the treatment of concrete problems, with comparisons to exact and / 
or independent full numerical solutions) of smoothed Wigner equations have already been 
examined in [2l[T], with very encouraging first results. 

The smoothed dynamics can be expressed as "convolution-deconvolution sandwiches" L = 

<I>L<I>"^ : ^w 1-^ ^Lw. (4) 

So the first question we treat in this paper is the explicit computation / representation of these 
"convolution-deconvolution sandwiches" , in ways amenable to analysis and computation. This 
comes together with the need for a basic framework, since the new operators are of infinite 
order in general. 

The core result can be outlined as follows (see Section [2] for the notations and conventions 
we use for the Fourier transform and the Weyl pseudodifferential calculus): 

Theorem 1.1 (Smoothed calculus) Let f{x) G >S(M"), L{x,k) G 5'(M^") and L be the 
operator with total Weyl symbol L{x,k), i.e. Lf(x) = f j^^^^e^'^^^^''^~'^^L{^^,k)f{y)dkdy . 
Denote moreover by $ the smoothing operator 

$ : f{x) ^ T-^\^ [e-f '^^^V.^x [/(a)]] , (5) 

^written here in a symbolic form, assuming that the operator L contains all initial/boundary conditions etc. 



(where T is the Fourier transform) and w = $/. Then the operator L = $L$ ^ can be 
expressed as 

L.(.) = 2" / e-<"-""=)—»(*-.)L,(2„. *)*<;*, (6) 

where 

Li{u,k) = J^^^u[Lix,k)]. (7) 

Another, equivalent formulation is 

Lw{x)= I e^'''^''-y^^~^^''''''L(^^,k]wiy-ia^k)dydk. (8) 



Naturally, before we can prove Theorem 11.11 we have to show that equations ([6]) and ([8]) 
make sense. This is achieved by showing that smoothed functions w{x) have properties very 
closely resembling those of Gelfand-Shilov functions of type S^^^ [T2], as we see in more detail 



in the body of the paper (section [33]) . 

A natural question to ask is "what is the Weyl symbol of a convolution-deconvolution 
sandwich $L$^^?" To motivate the answer, let us first consider the free Schrodinger evolution 



^^(x):=(e2 iIjq){x) = -== j e 2>* ^o{y)dy. (9) 

It is well known that the Wigner transform of a wavefunction ip is the Weyl symbol of its 
orthogonal projector, 

wmx,k) = aweyim{ij\). 

It is also well known that the free-space Schrodinger evolution ^ is pushed on the Wigner 
function level by 

W[ilJt]{x, k) = W[i]jq]{x + 27rtA;, fc), (10) 

which leads to 

W[tl)o]{x + 2'Ktk, k) = W[^t]{x, k) = e^^ iy[^o](x, A;) e"^^. (11) 

In fact this extends immediately to any Weyl operator composed by free Schrodinger evolution, 
i.e. in general 

c^weyi fe^^Le"^^j {x,k) = aweyiiL){x - 27rtk,k). (12) 

2 . 2 

By noting that $ = e^'^ which is formally equal to e*^'^ we can expect that, in the case 
of analytic symbol, we should have 

i (7 J> 
aweyi {^L^'^) {x, k) = aweyi (L) {x + ^-, A;). (13) 

This is in fact the case, as we prove in Section 13.11 



Theorem 1.2 (Smoothed calculus: Case of analyticity in x, Weyl symbol) LetL{x,k) G 
5'(R^"). Moreover denote Li{u,k) = JF^^^ [L(x, A;)] and assume that 3M > such that 
supp Li{u,k) C [—MjM]"^ X M", i.e. Li{u,k) has compact support in u. (In particular it 
follows then that for each k G M", L{x,k) is an entire analytic function of x). Then the Weyl 
symbol of $ L(a;, dx) $~^ is 

L{x,k) = L{x+'-^,k). (14) 

In the case that the symbol is analytic in k instead of x, we can also have a simplified version 
of Theorem 11.11 namely 

Theorem 1.3 (Smoothed calculus: Case of analyticity in A;) Let L{x,k) E S'(M?"). In 
addition assume that L{x, k) is a continuous function of (x, A;), and Vx G M" Lx{k) = L{x, k) 
is (the restriction to the real numbers of) an entire- analytic function, and moreover \fx G 
W^ G{k,y) = L i^^, k + ^^^~y' j g 5'(M'^). For example differential operators, L{x,k) = 

N 

^ Am{x)k"^, fall in this category. 

m=0 

Then 

<^L<!?-^wix) = [ F{x,k)e-^'"'''''wix-ia^k)dk, (15) 



where 



F{x, k)= [ e-2-(--^)^L(^, k + '^^^-yl)dy. (16) 



a2 



Finally it is important to emphasize that the explicit computation of the "sandwich" does 
provide us with important partial cancellations: applying $L$~^ as three distinct operators 
passes from a deconvolution, i.e. a Fourier multiplier of growth e'"'^ . Applying the result of 
Theorem 11.11 passes from imaginary translations, i.e. from multipliers of growth e'"^, but no 
deconvolutions. 

1.2 Homogenization in terms of the smoothed Wigner transform 

In this Subsection we briefiy outline the idea of SWT-based homogenization and its motivation, 
and state the main results in that direction. 

Consider a problem of the form 

ut + L{x,dx)u = 0, ,^^. 

u{x, 0) = uo{x), 



for a wavefunction u{x, t) : M"+^ -^ C^, where L{x, dx) is a. d x d matrix of pseudodifferential 
operators with matrix- valued Weyl symbol L{x, k). We will assume that equation (fTTl) describes 
the propagation of waves (e.g. Schrodinger, acoustics, Maxwell's equations etc). 
The physical observables of the wavefunction are the scalar functions of time 

M(t)= [ tF{y,t)Mu{y,t)dy, (18) 



corresponding to operators M = M{x, dx) from an appropriate class (e.g. with polynomials 
or Schwartz test functions as their Weyl symbols). Physically, bilinear observables describe 
e.g. energy and energy flux in many cases, including the examples mentioned earliecl. Auxil- 
iary quantities of interest in this context are the physical- and Fourier-space densities for the 
observables, 

M(x,t) = u{x,tfMu{x,t), 



.^{k,t) = u{k,tfMu{k,t). ^^^^ 



They are called densities because 



M{x,t)dx= / m:f{k,t)dk = M{t) (20) 

That is, equation flT7|) is seen here more as a book-keeping mechanism; the object of interest 
is not the point values of the wavefunction u{x, t), but a collection of observables and, to some 
extent, their densitieqj. This is often a satisfactory framework, most notably in quantum 
mechanics. 

The SWT-based homogenization approach consists in simplifying the book-keeping problem 
(fT7|) while keeping track exactly of the observables. The name hints towards the fact that the 
SWT results from the well-known Wigner transform, after the latter is convolved with an 
appropriate kernel. This is a very natural idea, and several variants have appeared in many 
contexts, most notably time-frequency analysis [TJ [16] and semiclassical limits (see the more 
detailed discussion after Theorem 11.51) . The SWT is defined as the sesquilinear transform 

W:f,g^W,,[f,9]{x,k) = 

We will work mostly with its quadratic (and time-dependent) version 

W[u]{x,k,t) = W[u{-,t),u{-,t)]{x,k) : r2"+i ^ C^x'^. 



^In most linear problems the natural energy functional is quadratic in the wavefunction; this can be under- 
stood better e.g. in the context of variational formulations. 

•^Another way to understand this point is that the problem (J17p describes the microscopic dynamics, but we 
may only be interested in a macroscopic view of the problem. 
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Despite the obvious increase in dimensionality introduced by the SWT, it can be used for 
compression, because it doesn't exhibit oscillations, in contrast to the wavefunction u. This is 
achieved by an appropriate smoothing in phase-space, controlled by the parameters ax, crj^. So 
the basic idea is switching an oscillatory wavefunction for a smooth phase-space density which 
lives on a twice-dimensional space. 

The simplest paradigm for this homogenization approach consists in two steps: 

• the derivation of exact equations for the evolution in time of iy['u](x, k,t), and 

• the derivation of a "smoothed trace formula" , expressing the bilinear observables directly 
in terms of Vr[M](x, k,t). 

(One more step is necessary for the treatment of systems, namely decomposing iy[u](x, k) on 
an appropriate matrix basis). 

The derivation of both the equations for the evolution of the SWT and the smoothed 
trace formula follows along the same lines as Theorem II. ![ and rests on the following core 
computation: 

Theorem 1.4 (Smoothed Wigner calculus) Let f{x),g{x) G S{W), L{x,k) G ^'(R^") 
and L be the operator with L{x , k) as its Weyl symbol. Moreover, denote w{x,k) = W[f,g]{x,k). 
Then 

W[Lf,g]ix,k)=jCW[f,g]ix,k) (22) 

where 

Cw{x, k) = 

^ 22n J ^2ni[Six-K+^a^,X)+nk+X+ialK)+xX+kK]^^g^ 2T)w{X -S,K- T)dSdTdXdK. 
X,K,S,T 

(23) 
or, equivalently 

Cw{x, k) = 

= J L{S,T)e'-'^'^'-^'^-^(-'^''+-'^^')w{x + ^^,k-^^)dSdT. (24) 



Like before, we can compute the Weyl symbol under appropriate assumptions: 

Theorem 1.5 (Weyl symbols for the smoothed Wigner calculus) Consider f(x),g{x) G 
5(M"), and L{x, k) to be the Fourier transform of a compactly supported tempered distribution. 
(In particular it follows that it is the restriction to the real numbers of an entire analytic func- 
tion). Then the operator C, defined in equation ^W) . has Weyl symbol 

C{x, k, X,K) = l(x- ^^^, k + ^^±f^) (25) 



^Calibrating the smoothing is pretty well understood, but it isn't central here. A rule of thumb is that a'^ 
must be comparable to the wavelengths of u{x, t), and af to the wavelengths of u{k, t); see also Section l4Jl and 
Section [5] for problems in the semiclassical scaling. 



in the sense that 

£.w{x, k)= J e^7rnx(x-a)+K{k^b)]l (^x±a^ k±b ^ ^^ K) w{a, b)dadbdXdK. ^26) 



aAX 



Observe how simple and intuitive is the passage to phase-space in terms of the Weyl symbols: 
Theorem 11.51 can be automatically guessed (and proved, if its assumptions hold) from Theorem 

The concept that certain bilinear functionals, and not the point values of the wavefunc- 
tion, carry the "important information" (the "physical observables" ) , originates in quantum 
mechanics (and has found applications in other contexts as well). Introduced in 1932 [23], the 
Wigner transform (WT) appeared in the 90's as an important tool for homogenization of wave 
propagation. The concept of semiclassical measures was extensively studied; see for example 
[T7t [Tot [HI m [25] ; adaptations to the case of Schrodinger operators with periodic coefficients and 
applications of the method to vector problems were carried out [U [201 E]) ^^d applications to 
stochastic problem were also studied, see e.g. [22]. 

In many of the works mentioned above the idea of smoothed Wigner transforms (often under 
the name "Husimi functions") appears as a technical device, e.g. for proving the positivity of 
the Wigner measure or for using interpolation estimates from classical kinetic theory, e.g. [19]. 
The link to coherent states (i.e. abstract wavelet transforms) has also been pointed out, e.g. 
in [T7]; however "a theory of smoothed Wigner transforms" has not been tackled, essentially 
because of the problem of dealing with "convolution-deconvolution sandwiches" and formulating 
explicit smoothed Wigner equations. 

This problem is solved here, and the formulation of exact smoothed Wigner equations for 
a broad class of problems, together with a smoothed trace formula for the recovery of the 
observables, is carried out in Section 14. 2[ As concrete examples, we work out the smoothed 
Wigner equations for the linear Schrodinger equation, the cubic non-linear Schodinger equation, 
and the Hartree equation. 

Virtually all the existing work with WTs is in the semiclassical regime, therefore it is 
appropriate that we look at a semiclassical application. We do so in Section [3 more specifically, 
we formulate and prove the following 

Theorem 1.6 (Semiclassical finite-order approximations to the smoothed Wigner calculus) 

Let N &N, V{x) : M" -> M. Consider a "semiclassical family of wavefunctions" {f^} C 5(]R") 
for which 3Mo > 0, £o ^ (0, 1) such that 

||r|U2(K.)^Mo V£G(0,£o). (27) 



According to Theorem\1.4 



/2 
e^^'s-''f-'s'V(S)w'ix + '-^S, k - -S)dS. (28 



Under appropriate assumptions for the potential V(x) |f|, and for a^ ^ 2, this expression can be 
approximated by differential operators 

W'[Vr,f%x,k) = E ^""Pm W'[f']{x,k) + r,{x,k), (29) 

m=0 



where Pm are homogeneous differential operators of order m, with coefficients depending on the 

n I /|2 

_ / \n 2-k\x — x I 

m- order derivatives of the smoothed potential V (x) = ( /^ ) J e ^ V{x')dx', and 



x'ei 



||re||L-(R2n) = C)(e^-"). (30) 

It should be pointed out that the - somehow surprising - possibly divergent bound of 
equation fl30l) has to be compared with the sharp estimate ||M^^[/^]||L°o(R2n) = O (^), as we 
will see in the case of WKB ansatz; see the first of the remarks after Theorem 15.81 in Section [51 

This result can be used to construct PDE approximations to the smoothed Wigner equations 
corresponding to semiclassical problems, in analogy to the construction of asymptotic equations 
for the Wigner measure in the works mentioned earlier (e.g. [IZIIII]). Such approximations were 
first proposed in [21 [1] , where it was seen that they provide an efficient computational method 
for semiclassical problems, recovering e-dependent information that the Wigner measure cannot 
keep track of. 

1.3 Proofs 

Theorem 11.11 is, verbatim, a concatenation of Theorem 13. 3[ equation (15T|) in Sections 13.11 and 
Theorem 13. 191 equation ( 113ip in Section [3^ 

Theorems 11.21 and 11.31 are Theorems 13.41 and 13.201 in Sections 13.11 and 13.41 respectively. 

Theorems 11.41 and 11.51 are Theorems 14.51 and 14.61 respectively, both found in Section 14. 1[ 

Theorem 11.61 is Theorem 15.81 of Section [51 

1.4 Organization of the paper 

In Section [3l we develop the necessary prerequisites, prove the core result and look into several 
equivalent formulations. We look into an elementary formulation, with minimal prerequisites. 



^Of the form 

3Ci,Mi > : \V{k)\ sC Ci\k\-^^' VO < |fc| s^ 1, 

and, there is an appropriate M2 < ~3 (satisfying additional constraints) such that 

3C2 > : \V{k)\ ^ C2\k\'^'^ V|fc| > 1. 



in Section [3m As we mentioned earlier, establishing certain properties of smoothed functions is 
a necessary step; this is done in Section 13. 3[ A very helpful tool (although not an exact match) 
comes from Gelfand-Shilov spaces, briefly reviewed in Section 13.21 Other formulations of the 
smoothed calculus are formulated and proved in Section 13. 4[ Indeed, these more sophisticated 
formulations will prove particularly useful in the sequel. The references to the specifics of 
Wigner transforms are kept to a minimum throughout these Sections. 

In Section 14.11 the fundamental calculus for the SWT is developed, as an application of 
Section [3l and a general-purpose phase-space reformulation scheme is outlined in Section 14.21 
An implementation of the phase-space reformulation of wave problems, as outlined earlier, is 
presented in Section 14. 2[ The smoothed Wigner equations for virtually any linear system, and 
the smoothed trace formula are formulated making use of the smoothed Wigner calculus. An 
interesting point is that closed phase-space equations for the cubic non-linear Schrodinger and 
Hartree equations can be obtained, essentially with no additional work. In Section[5]we examine 
the semiclassical asymptotics for the smoothed Wigner calculus. This allows for a quantitative 
comparison to the respective WT-based results. 



2 Definitions and notations 



The Fourier transform is defined as 



/>) = T.^k \J{x)\ 



e-'^"'''f{x)dx. 



(31) 



x& 



Inversion is given by 



T,^. {T-\ [/(a)]] = T^l, \Ta^, [/(a)]] = /(x). 
An operator L is denoted L(x, dx) and said to have Weyl symbol L(x^ k) when 

y,keR" 

An equivalent expression which we also use is 



(32) 
(33) 

(34) 



Lfix) = T-^lx 



u 



/(A;--)e"""Li(M,A;)rfM 



(35) 
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where Li{u, k) = Tx^u [-^(a^, ^)]- 
The trace formula, 



L{x,-k)W{lg\{x,k)dxdk= I Lf{y)-g{y)dy, (36) 

x,fceiR» 



can be seen as an equivalent definition of the Wigner transform (WT), 

W[f, g] {x, k)= I e-'-^'yf (x + I) ^ (x - |) dy, (37) 

associating it with the Weyl calculus. 

Sometimes a scaled version of the Weyl calculus is used. This can be motivated e.g. from 
WKB functions, i.e. functions of the form f{x) = A{x)e^^^^\ The scaled Weyl calculus is 
defined consistently with the scaled WT through the trace formula, i.e. L{x,edx) is defined 
through 

j L{x,edx)f{y)-g{y)dy= j L{x,k)W'[f,g]{x,k)dxdk, (38) 

yeR" x,keR" 

where 

WV,9]{x,k)= I e"'-^'yf{x+'-^)-g(^x-'-^)dy. (39) 

yR" 

It is often the case that the Weyl symbol itself depends on e, L'^{x,k) |j. In that case 
L^{x,edx) is defined as 

j L'{x,k)W'[f,g]{x,k)dxdk= j L%x,edx)f{y)-g{y)dy. (40) 

rc,A;eIR" yeK" 

This can easily be seen to be equivalent to 

L%x,edx)f{a)= I e'-^'^'^-y^L%^y,ek)f{y)dydk. (41) 

For a function F{x, k) satisfying appropriate conditions (and for any tempered distribution), 
the operator F{x, ^) is defined, in terms of the Weyl calculus, as the operator with Weyl 
symbol F{x,k). It must be noted that usually Weyl symbol classes are taken to be more 
restricted than S'. Imposing some more assumptions will probably be necessary in certain 
contexts; however for our purposes the more general choice works well. 



^e.g. L''{x, k) = P{x, k) + eQ{x, k). 
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The order S of an operator L{x,dx) is defined as the smallest S > s.t. Va,/? G N U 
{0}3Ca,(i >0 : V/,m = l,...,n 



Qa+I3 

dx^dki 



■L{xM<CaA^ + \k\f-^. 



(42) 



Finite order PDOs are well defined on Sobolev spaces of the same order [15]. Typical 
examples of infinite order operators include deconvolutions and imaginary translations. 



Remark on notation: Please note that our conventions and notations for the Fourier trans- 
form and the Weyl calculus, clearly stated here, are used throughout the text without additional 
explanation. 



3 Smoothed calculus 



3.1 Explicit formulation of convolution-deconvolution sandwiches 



In this Subsection we will present the derivation of the elementary formulations of Theorem ll.il 
(i.e. equations ([6]) ), focusing on the mechanics of the derivation, as well as basic interpretation 
and application issues. The pseudodifferential and other operator-theoretic aspects are kept to 
a minimum here. 



Definition 3.1 (The smoothing operator) The operator $ is defined as in equation ^, 



I.e. 



$ : fix) ^ T,l e^^^'^'^'T^.^, \j{x')\ = ^ / 



-2-K 



(x-x'Y 



f{x')dx'. 



xeM" 



(43) 



Remark: The notation e 2 °" '^ is us 
generalization of Definition \3. 1\ would he 



$.„...,.„ : f{x) ^ T-,l, 



ed interchangeably with e '2-" I*'' , i.e. k^ = k-k. A natural 



e ^- ' 'J'.'^kifix')] 



2T 



-27rE 



/ / -,2 



(44) 



;?i '^' xeK" 



/ e '=1 ''t f(x')dx'. 



In fact we will use a a smoothing like that later, but most of the time it isn't worth the notational 
inconvenience - and all our results are generalized to the anisotropic case in a straightforward 
manner. 
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Operators like $ are very common, and are often called mollifiers. Observe that $ is 
translation invariant, i.e. a Fourier multiplier. Indeed, its action is very intuitively seen in 
the Fourier domain (it damps "the high wavenumbers" with a Gaussian weight). It is also 
straightforward to observe that it is one-to-one, with inverse 

$-1 : /(^) ^ jr~^l^ [ef '^^^V.^.x [/(x')]] . (45) 

Understanding the image by $ of the Schwartz test-functions iS, as well as other basic spaces 
will also be important. In particular, it is fair to say that smoothed functions are restrictions 
to a real space of entire-analytic functions. 

Observe moreover that $~^ is very hard (often impossible) to implement in practice, e.g. 
numerically. This is the basic reason why we want to compute explicitly the "sandwich" $L$~^, 
looking for some sort of mutual (partial) cancellation of $ and <l>~^. 

Let us start with a very simple observation: 

Lemma 3.2 (Smoothed polynomial calculus) Vtti ^^^i^ {!,... ,n} 

$^r*-' = (^^ + ^j ' (46) 

$a^$-i = a^. (47) 

Proof: The way to interpret and prove any expression of the form $L$~^ = L is by checking 
that V/(x) G 5(R") 

L$/ = $L/. (48) 

Indeed, to prove equation fj46l) for tti = 1 it suffices to check that 



(49) 



0-" 47r 

n (x~x^'^ 

= ^ J e-'-^^x'f{x')dx' = ^xf). 

In order to prove equation ( 1471) it is easier to work in the Fourier domain: 

J'.^x [d.^^f] = e-t'^'^'27rXi/(X) = J^^^x [^d,J] . (50) 

The generalization for m > 1 for either case is obvious. The proof is complete. 
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Theorem 3.3 (Elementary formulation of the smoothed calculus) Let f(x) E iS(]R"), 
L{x, k) G 5'(]R^"), L be the operator with Weyl symbol L{x, k) and 

w{x) = $/(x). 

Then 

$L$-V;(,t) = 2" /" e2"(*^-"''=)"+2-*^^u;(A; - u)Li{2u, k)dudk, (51) 

where Li{u, k) = Tx-^u [-^(2^, ^)]- 

Proof: Let us check first of all that the integral of equation (ISTl) indeed converges. To that 
end observe that the rhs of equation ( ISTl) is equal to 



2n J e2^(^^-'"''=)"+2^*'=^u;(fc-w)Li(2n,A;)durfA; = 

fc,«GM" 

= 2" / e2^(*^-'"''^)"+2^^'=^e-t'^'('^-")'/(^ - M)^i(2ti, k)dudk = 

_2^ J e2'^"('=+")-f-'('=+«)'/>-M)^i(2M,A;)(iMdA;= ^ ^ 

fc,neM" 

= J Fx{k,u)Li{2u,k)dudk. 
It suffices to show that 

because if that is true, then the integral exists (for each x G M") as a duality pairing between 
Fxik,u) G 5(M2") and Li (2m, fc) G ^'(M^"). 
Denote 

T : F{k,u)^-^ F{k-u,k + u). (53) 

T is essentially a rotation, and it is clear that T (iS(]R^")) C 5(M^"). Now observe that 

F^(A;,ix) =2"r (^/(A;)e~^"'"'+2""") e SiR^""), (54) 

since Vx G M" 

/(A;)e-f"'"'+2™" G 5(M2") (55) 

and therefore equation ( l52l) can be cast as a Schwartzian duality pairing (the independent 
variables in equation ( l55l) are A;, m; x plays the role of a parameter). 

Morally, the point here is that despite the real exponential term in equation (ISTI) . the in- 
tegral exists because we act on smoothed functions (which have Gaussian decay in the Fourier 
domain) . 
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Let us now prove equation (!5T]) . Starting from the Ihs one observes that (we use equation 
(135|) from Section [2] for the implementation of the Weyl calculus) 



2"$ 






^x' 



k~*x' 



nn r „2iTi[xl—x'l+kx'+x'u 



J f{k- u) e2™"Li (2m, k)du 

f /(A;-M)e2™'"Li(2ii,fc)dM 



■J h-^'rl 



.U& 



\cjH' 



f{k — u)Li{2u, k)dudkdx' dl 



,,k,x'. 



(56) 



2^^ / e^nix'(~l+k+u)^^> g27r««-f a2«2^(^ _ ^)^^(2li, k)dudkdl 

u,k,x',leR" 

= T' J 5(A; + u-/)e2"^^'--'^'''' 



'd/ /(/c — u)Li{2u, k)dudk 



2n J e2'^^^('^+")-f'^'('=+")'/(fc-ti)^i(2M,A;)(iMdfc. 



M.fc.ei 



But we have already seen the last expression above to be equal to the rhs of equation flSTl) . in 
equation fl52l) . The proof is complete. 



It is very easy to check the consistency of Theorem 13.31 with the Weyl calculus for cr = 0, 
and with Lemma 13. 2[ Indeed the comparison with Lemma 13.21 reveals nicely the "miraculous 
cancellation" that takes place for polynomials: if L(x, k) is a polynomial, then supp^ Li{u, k) = 
{0}, and it doesn't allow the real exponential term to give rise to an infinite order operator. 
For example, if L(x, k) = {2mx)'^, equation (ISTI) becomes 



$L$ ^w{x) 



Mix-a^k)u^ U _ '^\ 5M(^ _ 0)d' 



u e 



2mkx 



dk, 



(57) 



fc,«eiR" 



generating of course the same end result as Lemma 13. 2[ In general however contributions from 
u away from zero will give rise to operators not of finite order (e.g. for L{x, k) a Gaussian). 
Observe finally that our result is, under appropriate conditions, equivalent to 



$L$-^ = Lix + 



0-2 d^., dr. 



4:71 27ri' 



(58) 



Theorem 3.4 (Weyl symbols for the smoothed calculus) Let L(x,k) E 5'(]R^"). More- 
over denote Li{u,k) = J-'^^u [L{x , k)] , and assume that 3M > such that supp Li{u,k) C 
[—M,M]"' X M", i.e. Li{u,k) has compact support in u. (In particular it follows then that 
for each k G M", L(x, A;) is an entire analytic function of x). Then the Weyl symbol of 
$L(x, (9^)$"^ is 

L{x,k) = L{x+'^-^,k). (59) 
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Remark: Observe that the Weyl symbols produced this way need not be tempered distributions 
(take e.g. L{x,k) = e^^). This is related to the fact that the smoothed operators are of infinite 
order in general. Understanding better what this means, both in terms of analysis and applica- 
tions, is one of the objectives of this paper. In any case it hints towards the use of more general 
test-function / distribution theories as a natural direction. We believe that in that context the 
assumptions on L{x, k) for this Theorem can be probably relaxed. 

Proof: We will directly verify that the operator with the Weyl symbol of equation fl59l) 
(equivalently the operator of equation ( l58l) ) is the same as that of equation ( 15T|) : 

Lw{x) = J e2"^'=(^-s')L (^ + ^, k) w{y)dydk = 

= 2" / e^^^'^^^-^^J^-j^ [e2™(^+2'+^'^'^Ui(2ii, fc)j w{y)dydk = 

(60) 



2n J e'"'['=("-^)+"("+^+'"''^)]«;(2/)dyLi(2n,A;)t/MrfA: 



3/,fceM" 



2" / e-2-^?'(^-«)«;(?/)di/e'"'N+"("+*"''^)]Li(2M,A;)(iMdfc 



= 2" / wik - u) e2-h+"(-+^-'^')]^^(2u, k)dudk. 
3/,fceiR" 

The first step which needs more explanation is the equality 

2-"L h^ + '-^,k) = T-^^ L2™(a.+y+*a2fe)^^(2«, A;)] . (61) 

To see that, first of all recall that L(x, A;) G iS'(M^"), and therefore Li(2n, A;) = Tx^v [-^(^^i ^)] \v=2u 
is well defined. Now 

/ e2™-(-+?'+-''=)Li(2M,A;)dii= / f l?^!!^:i^±fi!^Li(2n, A;)dM = 



(62) 

2_„ g [(-+^+-^^)-^-]' l(o, k) = 2'-L h-^ + ^:^,k' 



/=0 

The series is a Taylor expansion of L{x, k) in the x variable for each k; it converges absolutely 
following our assumption. This also justifies the interchange of the order of summation and 
integration in equation ( l62l) through dominated convergence. (See also lemma lS.lSp . 

Another step in equation ( 160|) that needs justification is the interchange of the du and 
dy integrations, passing from the second to the third line. There it suffices to observe that 
g-27ro- uk ^g^^ ]^g replaced by e~^''°" ^'^ X[-2M,2M]"{u), in which case the result follows by the 
standard tempered distribution calculus. 

The proof is complete. 
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As we coinmented briefly earlier, it seems reasonable that L{x + ^^^Y^,k) can be defined 
precisely, in an appropriate (weak) sense for any L{x,k) G 5'(]R^"'). Gelfand-Shilov spaces of 
test functions and their duals (ultra-distributions) seem like an appropriate framework in that 
direction. 

Before we get to that question, however, there are some more basic results about smoothed 
functions that we need; these too can be easily formulated as modifications of standard results 
for Gelfand-Shilov test-functions. In the following Section we go briefly over the basic facts 
concerning Gelfand-Shilov test-functions, to set the stage for our results concerning additional 
formulations of the smoothed calculus. 

3.2 A natural framework for the smoothed calculus: Gelfand-Shilov 
spaces 

Gelfand-Shilov spaces (and ultra-distributions) is a well-defined topic which has been attracting 
increasing attention recently, with several dedicated monographs and papers. It is clear that 
a full presentation of them is completely outside the scope of this work. What is however 
necessary, is to outline some basic facts and also some motivation as to why these - somewhat 
unusual - spaces are well suited for the study of our smoothed calculus. Here we will focus 
on presenting the necessary background; in the next Subsection we will focus on its use and 
application. 

Usual test-functions, i.e. Schwartz test functions, are defined by 

f{x) e SiW") ^ sup \xPdlf{x)\ < oo Vp,g G (NU {0})", (63) 

i.e., morally, by decay of any order of derivatives faster than any power of -ttti [12] • It is 
important to note that this condition is symmetric with respect to the Fourier transform, i.e. 

J^(5(M")) = 5(M"). (64) 

This allows for a theory of distributions "of some polynomial order of growth" in space and 
Fourier domain. In addition, the Weyl calculus allows the transfer of function-space and 
function-theory results to operators, providing very strong tools for many problems - in par- 
ticular problems involving differential operators. However, there are two classes of operators, 
very important in this study, that are not contained in this framework: imaginary translations 
(i.e. extensions to the complex plane) 

r,. : f{x) ^ fix + za) = T-^l^ p^'^V>)] , (65) 

and deconvolutions 

$-1 : /(x) ^ T,l^ [et'^^'=V>)] • (66) 
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So the basic motivation is very simple: we need to accommodate Fourier multipliers and / 
or Weyl symbols of faster than polynomial growth. (Of course things are actually more subtle 
than that in the end). This brings on the question of smaller test-function spaces - after all 
deconvolutions and imaginary translations are not well defined on all of iS(]R"). 

These spaces of "very nice functions" should then be examined with respect to (nontriviality 
first of all!), closedness under elementary operations, behaviour under the Fourier transform 
etc. This is contained in the Gelfand-Shilov theory of test-functions and distributions. 

At this point we will go over the Definition of Gelfand-Shilov test-functions and some basic 
facts that are relevant. All results are quoted from chapter IV of [12], and we use the same 
notation. The implications of these core facts for the problem at hand will be seen shortly. 
Many facts about them (after some adaptation) will be very helpful for working with smoothed 
functions. 

For simplicity we will work in the case n = 1, i.e. for functions of M. The generalization to 
higher-dimensional spaces is straightforward. 

Definition 3.5 (Spaces of type S) Let a,f] ^ 0. 

1. f{x) G Sa iff ^A,Cg depending on f s.t. 

\x''d^J{x)\ ^ CqA^k^'' VA;, g G N U {0} (67) 

2. f\x) G S^ iff 3B,Ck depending on f s.t. 

\x^&if{x) I ^ CuB'^q"^ Vfc, g G N U {0} (68) 

3. f{x) G 5f iff 3A,B,C depending on f s.t. 

\x^dlf{x) I ^ CA^B^k^'^q'^^ V/fc, g G N U {0} (69) 

Theorem 3.6 (Equivalent Definitions of spaces of type S) 1. f{x) G Sa iff^q G NU 
{0} 3a, Cq depending on f s.t. 

|a^/(x)|^C,e-"N". (70) 

The parameter a is related to A of the respective Definition 1 5*. 51 1 explicitly, and more 



specifically a = a{A) = — ^. 

eAct 

2. f{x) E S^ iff VA; G NU {0} 36, C^ depending on f s.t. f{x) can be (uniquely) extended 
to an entire function f{x + iy) satisfying the estimates 

\x'f{x + 2y)\^C',e''^y^^ . (71) 

The parameter b is related to B of the respective Definition \3.5[ 2 explicitly, and more 
specifically b = b{B) = ^{Be)^ . 
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3. f{x) G iSf iff 3a,b,C depending on f s.t. f{x) can be (uniquely) extended to an entire 
function f{x + iy) satisfying the estimates 



1 , ,^iv 



|/(x + 2i/)| ^Ce-'^l"l"+'l^l^. (72) 

The parameters a,b are explicitly related to the parameters A, B of Definition \3. 5i 3. 
Theorem 3.7 (Fourier transforms) 1. Fourier transforms of Sa 

^(5J = 5-. (73) 

2. Fourier transforms of S^ 

J'{S^)=Sp. (74) 

3. Fourier transforms of iSf 

T{Si)=Sl (75) 

Theorem 3.8 (Nontriviahty) The space S^ is nontrivial (i.e. contains a nonzero function) 
for all a > 0. The space S^ is nontrivial for all P > 0. The space S^ is nontrivial iff 

a + P ^ 1, a > 0, j3 > 0, or 

a = 0, P > 1, or 

13 = 0, a >1. 

Theorem 3.9 (Closedness under elementary operations) Let a,P ^ 0. Each of the 
spaces S^,Sa,S^ is closed under translation 



modulation 

dilation 

differentiation 

and multiplication by x 



f{x)^f{x + X), AgB 

fix) ^e'^^'-^ fix), AG 
fix) ^ /(Ax), A G M: 
fix) ^d^fix), 
fix) h^ xfix). 



Before we go on to formulate a more precise form of the smoothed calculus, we need to 
introduce a final family of spaces. 
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Definition 3.10 (Countably normed spaces of type S) Let a,P ^ 0, A,B>0. 

1. f{x) e Sa,A iff Vg G N U {0}, A> A, 30^^^ depending on f s.t. ^keNU {0} 

\x'dlf{x)\^C,,AA'k'", (76) 

or equivalently Vg G N U {0}, (5 > 0, 3Cg,5, 

\dlf{x)\<,C,,se~^<^^~'^\-\". (77) 

where a{A) = — ^. 

eAct 

2. f{x) G 5^'^ iff VA; G N U {0}, 13 > B ^C^^b depending on f s.t. Vg G N U {0} 



\X' 



'=a^/(x)|^C,,gSV^. (78) 



3. f{x) G 5^;^ iff \/A> A,B> B 3Ca,b depending on / s.i. Vg, fc G N U {0} 

|x'=a^/(2;)| ^ CA,BA''B''k'"'q''f^ (79) 

The new spaces are related to Sa,S^ ^S^ by 

SP = u 5^'«, (80) 

5q, = U iSqa, (81) 

A>0 

'^f = U 5f'f. (82) 

" A,B>Q "'"^ ^ ^ 

Moreover, they are ordered by 

Ai < A2 => 5„,Ai ^Sa,A2, (83) 

Bi < B2 => S^'^' C5^'«% (84) 

Si < B„ A,<A, => 5f;J C 5f;2. (85) 

Let us now summarize the generahzations of Theorems 13.71 13.81 13.91 for the new family of 
spaces here: 

Theorem 3.11 (Basic properties oi Sa,A,<S'^'^,SaA ) ^- Fourier transforms: 

.^(5„,a) = 5"'^. (86) 

.F (5^'^) =5^,5. (87) 

^ (<'f ) = 5;;^. (88) 
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2. Nontriviality: The spaces So,,a,<S'^'^ are nontrivial for all a,f3,A,B > 0. The space 
'^aA ^^ nontrivial iff 

a + (3 > I, a > 0, /? > 0, and A,B > 0; or 

a = 0, /3 > 1, A,5 > 0; or 

a > 1, /3 = 0, A,5 > 0; or 

a + /3 = 1, AB > 7(a, /9) /or some appropriate 7(a, /5) > 0. 

5*. Closedness under elementary operations: each of Sa,A,<S^'^ ,S^'j^^ is closed under 
translation, modulation, differentiation and multiplication by x i\. Dilations scale obvi- 
ously; if 

Dx : fix) ^ fiXx) 



then 



DxiS^^A) = S^^A, (89) 

DxiS^'"") = S^'^"", (90) 

DxiS'J) = 5f'f . (91) 



"'A 



3.3 Spaces of smoothed functions 

Remark: In this Section we will work in the general x G M" setup, in contrast to the previous 
Section, where we only examined n = 1. This is necessary, since the Wigner transform (which 
we want to apply our results to) doubles the number of independent variables, e.g. takes 1- 
dimensional problems to 2-dimensional ones. We go on to the more general n-dimensional case 
since it presents no essential difficulties. (Indeed Gelfand and Shilov also present briefly the 
n-dimensional generalization of their theory in Section 9 of chapter IV of IT^ . after a more 
detailed study of the 1-dimensional case). U 

Lemma 3.12 (The range of $, I) Consider a Schwartz test-function f{x) G iS(M"). Then 
its image under the smoothing operator belongs to a space of type S"^'^ , 



1 1 



f{x)eS{W) => ^f{x) eS^'^^% (92) 



^The operations are precisely defined in the statement of Theorem lS.QI 

^ For example, Definition 13.101 part 1, is generalized as follows: f{x) G Sa^(Ai,....A„) iff Vfc,(7 G 
(NU{0})", A =(Ai,...,A„), i,>^iV/ = 1,... ,71, 3C,_^ depending on/ s.t. 

1=1 

where of course x'^ = [xf'j , d^ = d'j.\. (Further generahzation with a — (ai,...,a„) is also possible, and 
straightforward, but it will not be of interest in this work). 
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or, more generally, 

f{x)eS => $.„...,.„/(x)g5^'^-(^'-'-), (93) 

Proof: We begin with the proof of equation (!92l) : we will prove that ^f{k) G Si i : then 

the result follows making use of Theorem 13.111 
Observe that, if we denote g{k) = $/(A;), 

g{k) = ^f{k) = f{k)e-^''''' =^ k^dlg{k) = e-^^'^''' J2 Pp,Uk) P (k) , (94) 

1=1 

for certain appropriate polynomials Pp^q-i{k) of degree at most p + q. f{x) G S implies that 

Wp,qeNU{0} 3Cp,, : \J2Pp,Uk)P\k)\^C,^g. (95) 

1=1 

Equations ( l94l) . ( l95l) together imply 

|Fa,^^(fc)|^C,,,e-i'^'l'=ll (96) 

Recall that equation ([77]) states g{x) G >Si _^ iff Vg G N U {0}, 6 > 0, 3Cg,5, 

2 ' ay/en 

|a|^(x)|^C,,,e-(t'^'-^)N^ (97) 

so equation ( 196|) is actually sirong'er than ^f{k) G 5i i . 

The operator $o-i,...,(t„ is defined in equation (jS]). Equation ( 193|1 follows in the same way as 
equation (192!) . 

The proof is complete. 

Observe also that, using equation (18^ . it follows that ^f{k) G Si i VA < a, and accord- 
ingly <l>/(x) G ^^'w^ VA < a. 

At this point it is clear that there is one more family of spaces we will use - and it is closely 
related to the Gelfand-Shilov test-functions: 

Definition 3.13 (Gi^,2, G'^'"') We will say that f{x) G Gi^^2 iff\fp,qe (N U {0})" 3Cp,g > 
s.t. 

\xPdU{x)\^C,,,e-^'^'\^\\ (98) 

and f{x) G G-2^^' iff 

f{k) G Gi^.. (99) 



2'" 
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Remark: In the anisotropic case a = (cti, ..., cr„) (i.e. the smoothing operator defined in 
equation (^^ ), we will use the same notation, namely: 
fix) G Gi,,. ^jfyp,qe (NU {0})" 3C,,, > s.t. 

|xPa^/(x)|^C,,,e"^£'^''', (100) 

and fix) G Gi'^' iff 

fik) G G.,,.. (101) 

We will not comment on the anisotropic generalization explicitly from now on, since it is 

n 

straightforward. It lies basically in the recasting of a'^x^ as ^ crfxj in the decay conditions. 

1=1 

The letter G is chosen to emphasize that these spaces are intimately tied to the Gaussian 
smoothing we use. 
Now it follows that 

Lemma 3.14 (The range of $, II) 

fix) G 5 ^ ^fix) G G^''^' ^ $/(A;) G Gi,,2 (102) 

and moreover 



2' 



VAi > a, < Aa < CT : ^^'^iv^ C G^''^ C 5^ 



(103) 



G2'- c 5^ 



where of course j- = i y—, ..., y^ j and similarly for the other indices. 



Proof: Equation (I102p essentially follows from the proof of Lemma [3 .121 Observe that equation 
dnHD (which holds for any /(x) G G^'^^) implies that if fix) G Gi ,^2 =^ e^^^fix) G S, and 

therefore, by a Fourier transform, if fix) G G^''^ then $^^/(x) G S. This shows $(5) ^ G^'"^ ; 
$(iS) C Gz'*^ is straightforward (in other words equations (!96l) and (198!) are the same). 
We will prove equation (11031) in stages; first of all we will show that 

VAi >(T : S^'^^ C Gi'^'. (104) 

Observe that (as we saw earlier in Definition 13.101 and equation ([9] 



fix) G 52'Mv^ ^ 
Vg G (N U {0})" , 5 > 0, 3G,,5 : |92/>)| ^ G,,5e-(i"?-^)^' ^ (105) 

Vp,gG(NU{0})",<5>0, 3G,,5: {x^d^ix)] ^ Gp,,,se-^^^'"-'^^\ 
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where Cp,g,5 can be chosen not larger than Cp^q^& ^ C s_ ■ sup 



\xr e ^2 



Therefore 



/(x)g5^'Mv^, Ai>(T => .^Qg^ 

Vp, g G (N U {0})" 3C;, : Ix^d^Jix)] ^ C^^e-f '^'-^ 
for C' ^ C a2-<t2 lI, and therefore, for Ai > a, S^'^^^^^ C G^'" . On the other hand 

p,q,-^-^ 

,— 1 7r2j2 12 — 1 

g'i(x) = ^j^_^^. e~2°" G Ga''^ \5^'^iv^, and now equation (11041) follows. 
That 

VA2 ^ a : G'5''^' C 5^'^^ (107) 

was shown in Lemma [3.121 (it was shown there for A2 = o", for A2 < cr it follows making use of 
equation (IMI)). Also, VA2 < cr : G2''^ C ^a'Aa^i^ jg obvious. 
The proof is complete lj. 

Lemma [3. 141 highlights the close relation of smoothed functions with Gelfand-Shilov spaces. 
This relation allows us to use (after some adaptation) a lot of existing theory - most notably 

1 2 

extensions in the complex plane. A very useful result has to do with the rate of growth of Ga'*^ 
functions on the complex plane: since $(5) C S'i'^ C 5^, estimates like the one in equation 
( !7T|) hold automatically. However, as we saw smoothed functions satisfy somewhat stronger 
conditions, and accordingly we can construct somewhat stronger estimates (which turn out to 
be necessary in the sequel). 

First of all however, we need a technical Lemma: 

Lemma 3.15 (Fourier-domain representation of entire functions) Let f{x) G 52(M'^). 
We know that then f{x) can be extended to an entire function on the complex domain, f{x + iy). 
In addition, the following representation is valid 

f(x + iy)= [ e^^'^-^''+'y^f{k)dk. (108) 



feeK" 



Remark: As we just saw (lemma \3.1^ smoothed test-functions belong inS^ ; therefore Lemma 
\3.15\ applies to them. 



9 



that is C' ^ Cp^qj where Cp^q^s are the constants of the same notation in equation (jlOSp for 6 = -k ^ „ 
^°It is not clear at this point whether 5^'°v^ is strictly larger than G's''^ , or if the two spaces are equal. 
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Proof: Denote z = x + iy; the proof lies with the computation 

OO 7 

/ e^^"'-^''+'y^f{k)dk= f ^ ''"^•if+^^)] f{k)dk 

fceK" 1=0 



OO 



I 
J2 ks{xs+iys) 



(109) 



= / E(2^0'^^^^^^ ^-fXk)dk = 

fceM" 1=0 

= / Ei'^mY'^t^— J, mdk = 

OO I OO 

= E E ^^'^ / k^fik)dk = E E f^5^/(o) = fiz). 

1=0 \r\=l fceK" «=0 \r\=l 

We have used the standard multi-index notation in the computations lj ■ 

In equation (11091) above it is first of all seen that the bulky expressions coming from the 
high-dimensional character of the problem can be nicely summarized as 



^2^ik-{x+iy)f(y^^j^_ f \^ 2mk\x + ly) 



f{k)dk= / Y. ^""'^Y'''\ fik)dk (110) 

«e{Nu{o})" 



etc. 

The condition we have to check to justify the interchange of summation and integration (in 
the more compact notation) is 



\/x,yeR: Yl *^V^* I \{2mk)' f{k)\dk<oo. (Ill) 



When we first expand in a series, in the first line of equation (11091) . it is the Taylor series of 
the exponential; when we summed analytically the series in the last line, this was the Taylor 
expansion of an entire-analytic function. That is, all the series converge, and the commutation 
of the series and the integral follows from the dominated convergence Theorem. 

Remark: Lemma lS . 1 5\ essentially is already justified for any entire-analytic function. However, 
since it is not too long, and to provide more insight, we will give here a more detailed proof of 
the result. In particular this highlights how the estimates on the entire function are inherited 
in this series - something we will return to. 



^^ re (NU{0})" is a multi-index; |r| =ri + ... + r„, r! = ri! • ... • r„!, k'' = k^ ■ ■■■■k';^", z"- = {x"^ +iy\^\ 

25 



Observe that f{x) G 1S2 ^ /(A;) g Si, (according to Theorem 13.71 part 2), and therefore 
(according to the high-dimensional version of Definition 13.51) . 3C,A = {Ai, ...,v4„) > : V/ G 
(NU{0})" 



{2Tnky f{k)\ ^ CAH^ = C n A';// 



;ii2) 



Now we have 



\x+iy\^ 



Ze(NU{0})" fcGM" 



«e(Nu{o})" 



/ \{2'Kik)^ f{k)\dk+ J \{2nikf f{k)\dk 



< 



^ Ciel^+^^l + E 



|a^+»;;l' 



ie(Nu{o})' 



/ \{2mkf^''>f{k)\j^^dk 



113) 



< 



^ Ciel^'+^^'l + C2 






I'^^S^ 



E ^^^('+2)(/ + 2)^, 
«e(Nu{o})" 



where in the last step we made use of equation (11121) . So now we only have to check the absolute 
convergence of the series 



E 



\x + iy\ 



-yl('+2)(/ + 2)^ = A^ J2 ^-^^^\ix + ^y)A\\ 



;il4) 



«6(NU{0})" Ze(NU{0})" 

or, equivalently, find the radius of convergence of the power series 

n 

E aiz^ = E n a/, 2;'= ; 
i&N" ie{nu{o})" *=i 



(115) 



1+2 



ls + 2 



ai 



n 

s=l 



ls\ 



to be infinite. 



Case 77, = 1: Observe that, making use of the Stirling approximation, it follows that 

^(v^)'(/ + 2)i(/ + l)i(/!)^, 



{1 + 2 



1+2 /(/ + 2)!e('+2) 



^/2n{l + 2) {2tt)- 



;ii6) 
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and therefore, using the Stirhng approximation once more, 



rii ;^-: - 1^- - ■ ■ ■ ^>', 



(117) 



^(( + 2)l(( + l)i(^). 



(27r)' 

It is now obvious that Vi?o > 3C = C{Rq) such that 

|a/|^C-/?o^- (118) 

But an estimate hke this imphes that the radius of convergence for the power series with coeffi- 
cients ai is at least Rq; therefore the radius of convergence is infinite, and the series of equation 
(11141) always converges. 

General case, n G N: This is a straightforward generalization of the previous computation. 
First of all, using the Stirling formula like earlier it follows that 

a.-n^^(/. + 2)l(/, + l)^(^)\ (119) 

and therefore Vi?i, .., -R„ > 3C = C{Ri, ..., Rn) such that 

|ad ^C- n R;^\ (120) 



It is a standard (and easy to show) Lemma that equation (I120p implies that the power series 
with coefficients a; (i.e. the series of equation (11151) ) converges whenever l^;^! < Rg. Since the 
i?s's can be chosen arbitrarily, it follows that the series of equation (I114p converges always in 
the multidimensional case as well. 
The proof is complete. 

Let us also remark that integrals like the one of equation (llOSp have also been studied under 
the name two-sided Laplace transforms 



Theorem 3.16 (Behaviour of smoothed functions on the complex plane) Let f{x) E 

1 2 

G'a''^ . As we saw earlier, it can be extended to an entire function on C"; moreover 

gix,y) = e-^y"fix + ty) E ^(M^"). (121) 

Proof: First of all observe the following elementary identity 

J',^, [i2mxrdlfix)] = i-dkf i-2mxr />). (122) 
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We will use the observation of equation (11221) and Lemma [3. 151 to prove equation fll25p . We 
begin from a slightly different point, i.e. 

{27rixf (27rz(x + iy)y d^^d'Jix + iy) = 

= {2iri{x + iy)f d'y f e'^^'''^''+'y\-dkY{-2Tiiky f{k)dk = 

fceK" 

= {27Ti{x + iy)y J e^^'^'~''+''y\-2'KkY{-dkY{-2'Kikyf{k)dk = 

fceM" (123) 

= / [dle^^'^^'^+'y^] {-2'Kky{-dkY{-2mkY f{k)dk = 



= {-ly J e^'''^^''+'y'^dl{-2'Kky{-dkY{-2mkyf{k)dk. 

It follows now, making use of the Definition 13.131 of G^'"" , that there is a constant Cp^^^,. such 
that 

\{2mxyi2mix + zyy''dld;fix + ty)\ ^ C^,,,,,, / e~'-'^-y~^'^''->^dy = C'^^^^^J^) e'^^y^ 

(124) 
The commutation of the dk integral and the dk derivative in equation (11231) follow from the 

1 2 

dominated convergence Theorem and the bounds for G^''^ functions. 

It is obvious how the following inequality follows from equation (11241) : 
\f p,q,r,s e {NU {0})" ^Cp^g,r,s such that 

I {27nxy {27nyy dldlf{x + ly) \ ^ C^,,,,,, e^^" (125) 

In order to complete the proof, observe that for any differential operator with polynomial 
coefficients P, equation (I125p implies that there is a constant C(P) > such that 

|P/(x + iy)\^ C(P) e^'y\ (126) 

Observe in addition, that for each differential operator with polynomial coefficients P there 
exists a different differential operator with polynomial coefficients P such that 



P(e-^(7(x,y)) =e-'"i^f{g{x,y)). 



Moreover if P is of order s as a differential operator, and its coefficients are polynomials of 
degree up to t, then P will be still of order s as a differential operator, and its coefficients will 
be polynomials of degree up to s + t. 
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Setting P = {27iixf {2myy d^d'y we get 

\¥{e~^y^f{x + iy)\ I =e''^y^\¥{f{x + iy))\ ^ C(P) (127) 

The proof is complete. 

1 2 

Corollary 3.17 Let f{x) he a smoothed test- function, i.e. f{x) G G'a''^ . Take a fixed x G M"; 
then 

g{y) = e-'-^'y'f{x + ia'y) G S (128) 

as a function of y. Moreover, for any fixed y, 

g{x) = f{x + ty) G S (129) 

as a function of x. 

3.4 Equivalent formulations of the smoothed calculus 

First of all, let us make a remark concerning the Weyl symbols for the smoothed calculus: 

Lemma 3.18 (Imaginary translations of distributions) Let L{x) G S' , y E M". Then 

1 2 

L{x + iy) is a well defined functional on G^'"^ . 
Proof: Take f{x) G G^'^'. Now 



L{x + iy)f{x)dx = / L{x)f{x — iy)dx, (130) 

which is well defined since f{x — iy) G S, as we saw in Corollary 13. 171 The proof is complete. 

So, basically, the idea is that since we act on "very nice functions" we can have more 
operations on our distributions, which will be interpreted weakly. Observe that the point of 
Lemma [3. 181 has absolutely nothing to do with actually extending L{x) into the complex plane. 
Giving meaning to the Weyl symbols of convolution-deconvolution sandwiches is somewhat 
more interesting, since we don't take a fixed imaginary translation, but go over an imaginary 
axis. 

In any case, it seems tempting to ask whether L(x + ^^^-^,k) simply belongs to an ultra- 
distribution space (i.e. to the dual of some space of the type e.g. S^'^). It seems probable that 
smoothed operators can be cast in a satisfactory framework simply as operators with ultra- 
distributional Weyl symbols. 

Now we go on to the equivalent formulations of the smoothed calculus, making use of the 
properties of smoothed functions that we just proved. 
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Theorem 3.19 (Smoothed calculus) Let f{x) G S(U.), L{x) G ^'(IR^")^ L be the operator 
with Weyl symbol L{x, k) and w{x) = $/(x). Then 

$L$-S/;(x)= f e''^'^''-y'^^-^^'''^'L('^-^,k\w{y-ia^k)dydk. (131) 

Remarks: Before we go on to the proof, some comments should be made: 

• Equation I^131\) is well defined, through Theorem \3.1(A 

• This form makes clear what we gain by computing explicitly the convolution- deconvolution 
sandwich, as opposed to applying $L$~^ successively as three different operators: in 
order to implement the sandwich we need to compute / implement imaginary translations 
w{x) 1-^ w{x-\-iy), which correspond to a Fourier multiplier e'^'^^^ , but not deconvolutions, 
which correspond to a much stronger Fourier multiplier, e^" ^ . 

Proof: It suffices to show that equation (11311) defines the same operator as equation ( 15T1) . 

^m-^w{x) = T I e^Ai^-''^k)u+2mkx^Q^ _ u)L^(2u, k)dudk, (132) 

k,u£R" 

where always Li{u, k) = Tx~^u [-^(^^i ^)]- We start from the expression of equation (I13ip : 

J g2^i(x-j/)fc~2^a2fc2^ (^x±y^ ^^ ^^y _ ^^2-^Yydk = 



^ J e^'''^^^y^''-^'"'''''L{^,k)Tu^y e2--''=«w(M) dydk = 

y,k£R" '- -' 

^ J ^2^i{x-y)k-2^a^k'^+2^iuy+2na^kuj^ (^£±J^^ ^^ w{u)dudydk = 

u,y,k&R" 

^ J e-'^^^y^^-''^L{^,k)dye'^-'^''^+'^^''''^^''-^^w{u)dudydk= (133) 

= T' J e2'^^('=-")^Li(2(fc-^/),A;)e2'^*^'=+2'^'"'^("-^)u;(u)(iMdy(iA; = 

_ 2" J ^2n^[ux+xk+^a^uk]l^ ^2^,, A;) w{k - u)dudydk = 

u,y,keR" 

= ^L^^^w{x). 

One more equivalent formulation exists when the Weyl symbol L is a differential operator, 
which is somewhat simpler: 
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Theorem 3.20 (Smoothed calculus, a reformulation for differential operators) Consider 
the same assumptions for f , L as in Theorem \3.19[ above, and in addition let us suppose 
that L{x,k) is a continuous function of {x,k), and Vx G M" Lx{k) = L{x,k) is (the restric- 
tion to the real numbers of) an entire- analytic function, and moreover "^x G M" G{k,y) = 

L i ^^1 k + ^i^^ j G iS'(M"). For example differential operators. 



N 



L(x,A;) = 5^A„(x)A:'", 

m=0 

fall in this category. Denote also w = $/. Then 

<!?L<!?-^w(x) = I F{x,k)e~^'"'''''wix-ia^k)dk, 



(134) 



where 



F{x, k)= [ e-2-(-^)^-L(^±^, k + ^^^^ 



)dy. 



(135) 



(136) 



Proof: It is clear that under our assumptions the statement of the Theorem makes sense. 
For the proof, it suffices to make the change of variables 



k = k 



I , K^-y) 



a^ 



in equation (11311) : 

$L$-i«;(x) = / e2-^(^-2^)^-2--''='L (^, A;) w{y - ia^k)dydk = 

= J e^-(^"^)(^'+^^)-2-^^(^+^^)'L(^,fc + ^)dyt.(x-za2/t)rfA;'= (137) 
_ J ^-2.i{x-y)k^ fx^^ ^ ^ l{^^ dye-^^''^^^w{x - ia^k)dk. 

y,k&R^ ^ ' 

If we have a differential operator as in equation (11341) the last expression is equal to 



<I>L<I> ^w[x) = J 






-2Tva^k^ 



w{x — ia^k)dk. 



(138) 



The proof is complete. 
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4 Smoothed Wigner homogenization 



4.1 Smoothed Wigner calculus 

In this Subsection we will derive the smoothed Wigner calculus, which, as we briefly described 
in Section 11.21 allows for the derivation of smoothed Wigner equations and a smoothed trace 
formula. This work essentially follows the same lines as Theorem 13.31 and its proof, being 
somewhat more complicated due to the specifics of the Wigner calculus. 

Definition 4.1 (The Wigner transform) We define the Wigner transform (WT) as the 
sesquilinear transform 



J e-'-^'yf (x + I) ^ (x - I) dy. (139) 



W:f{x),g{x)^W[f,g]ix,k) 



The generalization to vectors is straightforward, i.e. if f{x),g{x) : R" ^ C^ then 

[W[f,9]i^,k\j = W[f„g,]{x,k) i,j e {1,...,4 

The WT is well defined and continuous as a bilinear mapping 

W : 5(M") X 5(M") ^ S{ 
W : L'^(W) X L^(W) ^ L2( 



[)2ra\ 
i)2n\ 



(140) 

(141) 
(142) 



The WT has a number of properties which allow the interpretation of its quadratic version 
Vr[w](x, /c) = W[u,u]{x,k) (often called the Wigner distribution of u to avoid confusion) as a 
"time-frequency energy quasi- density" . That is 



W[u]{x,k)dxdk (143) 

is somehow "proportional to the energy (L"^ norm density) corresponding to the to the wavenum- 
bers k at the locations x for [x, k) G A". Making precise this interpretation (and understanding 
its limitations) is a classic topic in time-frequency analysis [?iLl3j|, and there is no need to stay 
on it too long here. 

One of the first findings however, is that the WT exhibits so-called "interference terms" , 
i.e. fast oscillations in phase-space, which severely limit its numerical and intuitive use. The 
interference terms are due to the non-linearity of the transform; for example in certain many- 
component signals (such as finite sums of Gaussian wavepackets) the "bad terms" can be exactly 
isolated as the cross-terms, 



W 



E 9n 
m=l,...,M 

+2Re 



{x,k)= E W[g„,]{x,k) + 

m=l,...,M 



E T.W[gp,gq]{x,k) 

p=l,...,M q<p 



(144) 
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In most cases however, isolating explicitly the "bad part" is not possible; the term "auto- 
interference" is used to emphasize that. The oscillations in phase-space are in general at least 
as fast as the oscillations in u (i.e. comparable wavelengths), but can be arbitrarily faster PI. 
This makes absolutely necessary some step of regularization; indeed in most applications of the 
WT some additional regularization device is proposed, be it convolution with a smooth kernel 
(similar to what we do) [16j, an appropriate scaled limit (in which the oscillations vanish) 
[T7t [TT] . or the introduction of a stochastic averaging [22]. For a more complete discussion of 
the WT's interference terms and their interpretation, the interested reader can see [T^f [7]. 

Definition 4.2 (Smoothing in phase-space) Denote by $ the operator 
$ : w{x, k) ^ Tt,K^,,u [e"^['''''''''^^''Va,fe.x,i^ [«)(a, h)\ = 
= ^ / e ^^ ^^w{x',k')dx'dk'. 



xeif 



We use the same symbol as in Definition 13. ![ although technically it is a different opera- 
tor. Still, we will go on with this abuse of notation, because they are essentially very similar 
operators, and it is very easy to understand which one is used from the context: the one of 
Definition 13.11 acts on functions of x e M", while the one of Definition 14.21 acts on functions of 

{x, k) G M^". 

Definition 4.3 (The smoothed Wigner transform) The SWT is the sesquilinear trans- 
form 

W:f,g^W[f,g\{x,k) = 

= (£) / ^ ' ^^f{u+l)-g{u-l)dudy= (146) 

= ^WU\g]{x,k). 

The generalization for vectors is the same as for the WT. 
Moreover, it is well defined as a bilinear mapping 

W ■ SiW) X S{W) ^ G'^'('^-'^fe)(M2n) c 55'^^^'^\m2") c 5(M2"), (147) 

W ■ L'^iW) X L'^iW) ^ L2(R2"). (148) 

Of course we will be working a lot with the (quadratic) smoothed Wigner distribution 
W[u\{x,k). The parameters CTx^f^k control the length scales of the smoothing. The motiva- 
tion is to smooth out any oscillations at length-scales finer than the oscillations of u{x) itself 
{a1 is scaled with them) or those of u{k) {al is scaled with these). It should be mentioned here 
that if 

a^ak = 1 (149) 



^^Take f{x) = e ^ ^gix) — e '^^ "-^ . Then W[f + g]{x,k) has oscillations with wavelengths of order -, while 
the function f + g itself is not really oscillatory at all. 
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then 

W[u]{x,k) = (^j ' I I e-^^^^^-^--- uiy)dy\\ (150) 

i.e. W coincides with a spectrogram (also known as Husinii transform) with window g{y) = 

— o 

( ^ ) e ^' , and is therefore nonnegative. We will say that when axCTk = 1 we have crit- 
ical smoothing, while if a^cCTk < 1 the smoothing is sub-critical. Generally speaking, critical 
smoothing is pretty strong, and over-critical choices a^cTk > 1 are not interesting. So (jl,al 
are measured in units of typical wavelengths of u{x), u{k) respectively, and the strength of the 
smoothing is gauged by the number (Jxcyk G (0, 1]. This automatically puts some structure in 
the parameter space, which is found to be sufficient in many practical applications - although 
clearly there is room for more quantitative results in this respect. For more discussion and 
examples on the calibration of the smoothing see [21 [1] . 

Before we go on to the smoothed Wigner calculus, let us formulate, in our notation, the 
Wigner calculus: 

Lemma 4.4 (Wigner calculus) Let f{x),g{x) E 5(M"), L{x,k) E iS'(]R^") and L the pseu- 
dodifferential operator with L{x, k) as its Weyl symbol. Then 



where 



W[Lf,g]{x,k) = CW[f,g]{x,k) (151) 



Cw{x, k) = 

= J e2'^*K+HL(a, b)w (x + f , A; - |) dadb = 

(152) 



92n'r— 1 

^ ''X.K^x,k 



J ^2ni[s{x-K)+T{k+x)]^2S, 2T)w{X -S,K- T)dSdT 



or, equivalently, its Weyl symbol is 

C{x,k,X,K) = L{x-—,k + —) (153) 

We state Lemma S3] for completeness and motivation; it is a standard result, and the proof 
is also contained as a special case of Theorem 14. 5[ 

Theorem 4.5 (Smoothed Wigner calculus) Let f{x),g{x) E S{W), L{x,k) E ^'(M^") 
and L the operator with Weyl symbol L{x, k). Then 

W[Lf,g]{x,k) = CW[f,g]{x,k) (154) 
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where 

Cw{x, k) = 

^ 22n J ^2ni[S{.-K+^alX)+Tik+X+ialK)+.X+kK]l^^g^ 2T)w{X -S,K- T)dSdTdXdK, 
X,K,S,T 

(155) 
or 

Cw{x, k) = 

= J L{S,T)e''''^'^-'^'^-^(-'^''^-'^^')w{x + ^^,k-^^)dSdT. (156) 

Remark: First of all, let us remark that Theorem \4.5\ can be seen as an application of The- 
orem {37^ using the Wigner calculus W[Lf,g]{x,k) = CW[f,g]{x,k). In that connection, 
£ = ^C^^^. However, we will prove Theorem \4 . 5\ similarly to, but nevertheless independently 
from Theorem \3.3l one reason is that computations which are anyway necessary when working 
with SWTs will be carried out in the process. 

Proof: First, we will see that the operator C is well-defined on G^'^'^^''^k)(M?^) functions for 
each of the formulations of equations fll55p . fll56p . As we saw in Theorem 13. 16[ 



(157) 



9iiS,T) = 

= e^-^(3-+T''yH'^'^^'+4T-)uj{x + ^^, k - '-^) e 5(M2n) 

as a function of {S,T), and therefore equation (11561) makes sense. 

To see that equation (I155P is well defined we have to demonstrate that 

9{S,T) = 

^ J e2-[^(--^+-?^)+^('=+^+-^^)+-^+'=^]u;(X-5,K-T)rfX(ii^G5(M2n). (158) 

X,K 



But 



9{S,T) = 

= r ^2ni[Six-K+ialX)+Tik+X+ialK)+xX+kK]^l^j^ _ S,K - T)dXdK 



X,K 



2m[2Sx+2Tk]-2TT\crlS'^+crlT'^'\ fl ^Ql 



X,K 



^ ^2.i[2Sx+2Tk]-2.[.lS-+.lT^]^^^ ^ j, ^ ^^2^^ ^ _ ^ ^ ^^2j.y 

The last equality makes use of Lemma [3]T5l The end result is a Schwartz test-function according 
to Theorem 13. 16[ 
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We will show that equations (11551) . (11561) are equivalent. Indeed, the passage from equation 
(11551) to (11561) is essentially demonstrated in equation (I159p . Observe that we only do a change 
of variables and a Fourier transform, so the reverse course follows as well. 

So we checked that all the formulations in the statement make sense and are equivalent. 
Now we will finally show that they give the smoothed Wigner calculus, i.e. equation (11541) 
holds. Like earlier, the way to check it is by showing 

<l>Cv{x,k) =C<l>v{x,k), (160) 

where w{x, k) = ^v{x, k), i.e. v{x, k) = W[f,g]{x, k). The Ihs of equation (I160p is equal to 

^Cv{a, h) = <^L{x - Ife, A; + ^)v{a, h) = 



r\2n -T — 1 

^ -'X.K^x,k 



J ^2ni[s{x^K)+T{k+x)]l^2S, 2T)v{X -S,K- T)dSdT 



S,T 



(161) 



= 2^" r ^2-Ki[S{x-K)+T(k+X)+xX+kK-xA-kB+aA+bB\-^[a-lA^+crlB'^) 

L{2S, 2T)v{X -S,K- T)dSdTdXdKdxdkdAdB = 

= 2^" C^2m[-SK+TX+x{X+S-A)+k{K+T-B)+aA+bB]~^{alA^+alB^)^^^^ 

L{2S, 2T)v{X - S,K - T)dSdTdXdKdAdB = 

= 22" j-5{X + S- A)5{K + T- B)e^^'l-^^+^^+-^+^^yU-'^^'+-l^')dAdB 
L{2S, 2T)v{X - S,K- T)dSdTdXdK = 

= 2^" r ^'i^i[-SK+TX+a{X+S)+biK+T)]~^{a'^{X+Sf+al{K+T)^) 

L{2S, 2T)v{X -S,K - T)dSdTdXdK. 
The rhs of equation (11601) is equal to 
C,^v{a,h) = 

= 2^" C Q'^-^i[S{a-K+icrlX)+T(b+X+icrlK)+aX+bK]^2g 2T) 

^-^-{alix-s?+aliK-Tf)^^X - S,K- T)dXdKdSdT = 

^ 2^" C ^2m[S{a-K)+T(b+X)+aX+bK\j^t2S 2T) i^^"^) 

^-^-(.Ux+sf+aliK+T)%^X -S,K- T)dXdKdSdT = 

^ 22" r g2wi[~SK+TX+a(X+S)+biK+T)] 1(^23 2T) 

^-^-{.lix+s?+al{K+T?)^(^X -S,K- T)dXdKdSdT, 
which is the same as the last member of equation (11611) . 
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The proof is complete. 



Theorem 4.6 (Weyl symbols for the smoothed Wigner calculus) Consider f{x) , g{x) G 
5(M'^), and L{x, k) to be the Fourier transform of a compactly supported tempered distribution 
(in particular it is the restriction to the real numbers of) an entire analytic function. Then the 
operator C, defined in equation ^154^ , has Weyl symbol 

C{x, k, X,K) = l(x- ^^^^, k + ^±fi^) . (163) 



Proof: This proof follows along the exact same lines as the proof of theorem 13. 4[ 
Indeed, observe that 

L{x-'-^^^,k+'-^^^)u.ix,k) = 

^ J ^2.i[.X+kK-aX-bK]^ (^^ _ K-iaiX^ ^ ^ X+i^^ ^^^^ b)dadbdXdK = 



2-Ki 



^ ^ + a-K + i,ylx ^^ k + b+X + zcrlK ^^ ^^ 



= / e L ^ ^ i L{s,t)dsdtdSdT 

^2ni[xX+kK^aX-bK]^^^^ ^-^dadbdXdK = 
^ J ^2ni[x{x+^+ial§)+K{k-§+ial^)+S^+T^^a{X-§)-b{K-^)] 

L{S, T)w{a, b)dSdTdadhdXdK = (164) 

^ J ^2ni[x{x+^+ial§)+K{k-§+ial^)+S^+T^] 

L{S, T)w (X - f , ir - f ) dSdTXdK = 

= 2^" r p2TTi[X(x+T+ialS)+K{k-S+i(7lT)+xS+kT] ^2Q 2T) 

w{X -S,K - T)dXdKdSdT = 

= 2^" r 2Tri[Six-K+ia^X)+T{k+X+ialK)+xX+kK] j^/2Q 2T) 

w{X -S,K - T)dSdTdXdK, 

which is exactly the rhs of equation (11551) . 

We have seen the justification of imaginary translations through Fourier transforms for entire 
functions before. Here it is applied on a smoothed function function and therefore Lemma [3. 151 
applies. 

The other step that needs justification is the interchange of the dadb and dSdT integrations, 
passing from the third to the fourth line. Remember that we have assumed 3M > such 
that suppL{S,T) C [— M, M]^", therefore the real exponential terms e-'^'^x^s-ira^KT ^^^ ^^ 
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substituted by e~'^°"^"^'^"'^°"fc^"^X[-2Af,2M]2"('S', T) without changing anything. The result then 
follows by the standard tempered distribution calculus. 
The proof is complete. 



4.2 Coarse-scale dynamics in phase-space 

Theorem 14.51 allows us to carry out in a precise manner the basic steps of SWT homogenization 
outlined in Section 11.21 



Corollary 4.7 (Smoothed Wigner equations) Let L{x,k) E {S'{M?'^)) ^ , L = L{x,dx), 
Uo{x) e (5(R"))'^. Consider the IVP 



Ut{x, t) + Lu{u, t) = 0, 
u{x, 0) = uo{x). 



(165) 



Then the SWT of u, 
satisfies the IVP 



Wix, k, t) = W[u{-, t)]{x, k) (166) 



Wt{x, t) + 2n (cw{u, t)\ = 0, 

W{x,k,0) = W[uo]{x,k), 



(167) 



where C is defined in terms of L as in Theorem \4.5\ and Ti-^A) = ^ii^ denotes the Hermitian 
part of a matrix. 



The proof is obvious, and consists in the application of Theorem 14. 5[ and the observation 
that 



^ W{x, k, t) = W[ut, u] {x, k, t) + W[u, ut]{x, k, t) 



at 



-W[Lu,u\{x,k,t) — W[u,Lu]{x,k,t) = (168) 

-jC,W[u]{x,k,t) - (jCW[u]{x,k,ti * 



Corollary 4.8 (Smoothed trace formula) Let Mij{x,k) G iS'(]R^"'), i,j = l,...,d. Any 
quadratic observable of a wavefield Ui{x,t) G 5(]R"), i = l,...,d, corresponding to the oper- 
ator M = M {x, dx) (defined as in equation D^) ) can be directly expressed in terms of the 
(Hermitian-matrix-valued) smoothed Wigner distribution of u as 

M{t)= I tr(MW[u]{x,k,t))dxdk, (169) 
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where Aiij is defined in terms of Mij{x, k) in the same way as C is defined in terms of L{x, k) 
in Theorem \4.5[ Moreover, observables can be resolved over phase- space at coarse-scale, 

M{x,k,t) = tr(MW{x,k,t)) (170) 

consistently with their natural resolutions, 

M{x,k,t)dk = -^ / e ^ tr{if{x',t)Mu{x',t))dx'. (171) 

X J 



/O" r „2,r(fc-fc')2 

{x,k,t)dx = - — / e ^'i tr(iF{k',t)Mu{k\t)]dk'. (172) 



and 



k J 

The proof is obvious, since 

MW[u]{x,k,t) = W[Mu,u]{x,k,t), (173) 

and, Vw,t;e (5(M"))'^ 

/ tT(w[u,v]{x,k)\dxdk= / v'^{x)u{x)dx, (174) 



vr 



W[u,v]{x,k)dx = ^ / e ^fe tr{¥{k\t)u{k\t))dk\ 



(175) 






and similarly for the dk marginal. 

In particular, all observables corresponding to polynomial Weyl symbols (which typically 
include energy and energy flux) can be recovered from the SWT of the wavefunction in terms 
of finite-order operators. 

Corollaries 14.71 and 14.81 show how the smoothed Wigner calculus allows us to reformulate 
problems, originally formulated for "waves" (i.e. for an oscillating wavefunction on M"), to 
problems for "phase-space densities" (i.e. smooth / simple functions on M^"). Indeed, many 
well known paradigms fit in this general description, with semiclassical limits and Wigner 
measures being the most relevant from a technical point of view [31 El El El [13 Ell 



The introduction of a "fundamental length", controlled by cr^, a^ is a distinctively different 
feature from the Wigner measure approach; the concept of a fundamental length has also been 
discussed from a physical point of view as well. For more treatments that can also be described 
as "phase-space homogenization" in a wide sense - but not as closely related to what we do here 
from a technical standpoint - see also [22j and the relevant survey in the introduction therein. 
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A different problem (although similar in the sense that it involves infinite-order equations 
governing a smooth density) is treated in [T8] . 

Naturally, a concrete example is in place here: 

Example 4.9 (Schrodinger equation) Consider a wavef unction satisfying the Schrodinger 
equation, 



% (x, t) - \Au (x, t) + iV (x) u (x, t) = 0, 



176) 



u{x, 0) = Uo{x). 
Then its SWT W{x,k,t) = W[u{-,t)]{x,k) satisfies the equation 

|iy (x, k, t) + {2nk ■ V. + 4 V. ■ Vfc) W (x, fc, t) + 

^/e V My(s)iy(x+^,fc-f)ds , n77) 

iy(x,fc,0) = iyK](x,A;). 

The algebra of the smoothed calculus yields a very fortunate "accident" : we can have exact 
coarse-scale reformulations for certain nonlinear equations with no additional work: 

Example 4.10 (Cubic non-linear Schrodinger equation) Consider a wavefunction satis- 
fying the cubic NLS equation, 

^M (x, t) - I Aix (x, t) + i {Vi (x) + p\u{x, t)f)u (x, t) = 0, .^^g. 

u{x, 0) = Uq{x). 

Then its SWT W{x,k,t) = W[u{-,t)]{x,k) satisfies the equation 

§-^W{x, k, t) + (2'Kk ■ V. + f V. ■ Vk) W (x, k, t) + 
+2Re[i J e2---F(s,t)W^(x + ^,A;-f)ds) , z^^gN 

W{x,k,0) = W[uo]{x,k), 

F{s,t)=J^,^s[F{x,t)], 

F{x,t) = ^ J e -^ Vi{x')dx' + p J W{x,k,t)dk. 



where 



(180) 



x'eM" fceK" 
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Example 4.11 (Hartree equation (to smoothed Vlasov)) Consider a wavefunction sat- 
isfying the Hartree equation, 



:i8r 



f^u{x,t) - ^Au{x,t) + ilVi{x)+ f Vo{x-x')\u{x',t)\'^dx' \u{x,t) = 0, 



x'eii 



u{x, 0) = uo{x). 
Then its SWT W{x, k, t) = W[u{-, t)]{x, k) satisfies the equation 



at 



^Wix,k,t)+{2nk-W, + ^V^-Vk]Wix,k,t) + 



where 



W{x,k,0) = W[uo]{x,k), 
F{s,t)=J^,^s[F{x,t)], 



(182) 



F{x,t) = ^ J e -^ Vi{x')dx' + P J Vo{x-x')W{x',k,t)dkdx'. 



183) 



The derivation for either nonlinear equation follows by observing simply that the potential 
appears in equation f ll77p not just as V{s), but as 



e-f-^^V(s) = J-,^, V{x) 



•J X- 



vr 



a'' 



'V{x')dx' 



x'& 



i.e. instead of the original potential V{x), it suffices to know the smoothed potential, V{x). 
This, coupled with the marginals property of the SWT - equation fll7ip for M = I - makes it 
possible to have closed smoothed Wigner equations in this case. 

As we mentioned earlier, the closest relative to this approach (and an important motivation 
for it) is the WT / Wigner measure based semi classical limits technique. In the next Section 
we study an application of the SWT to semiclassical problems. 



5 The semiclassical regime 

Let us start with a few words of motivation. We will work in an asymptotic regime, scaled with 
a parameter < e << 1. The intuitive meaning of the small parameter e is that we work with 
signals / functions which exhibit very fast oscillations, e.g. WKB functions 



r{x) = A{x)e 



^^S{x) 



;i84) 
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Under certain conditions (e.g. the envelope A{x) is itself "smooth", "slowly varying") it can 
be said that the function of equation (11841) has amplitude n{x) ~ |y4(x)p and "instantaneous 
frequency" / "local wavenumber" k{x) ~ VS{x). Indeed abstractions like these - and making 
them precise - are at the heart of the (motivation for the) WT and time-frequency analysis. It 
must be clear already why it is natural that these questions are formulated in an asymptotic 
regime e << 1. Fi 

Definition 5.1 (Semiclassical scaling of the WT) The semiclassically scaled WT is de- 
fined as 

W : f,g^W'[f,g]{x,k) = 
= J e~'-^'yf{x+f)-g{x-f)dy = ^W[f,g]i!l). (185) 

in agreement to (T^ [721 IMS- 

Definition 5.2 (Semiclassical scaling of the SWT) The semiclassically scaled SWT is de- 
fined as 

W^:f,g^W%f,g]{x,k) = 

2 2 o 

^) / e-'""'''' ^ -^^f{u + f)-g{u-f)dudy= (186) 

= i^-^Vi^.v^.. W^[/, 9] {x, f ) = $v/i..,vi., w^l/, g] (x, k) . 
Obviously the choice of the scaling of the smoothing, 

0"a;, CTfc t-^ a/£0-^, ^/sak-, (187) 

is to some extent arbitrary; in what follows we hope to show that it is a natural choice, at least 
for some problems, lj 

A central object in semiclassical problems is the well studied Wigner measure (WM). We 
recall (adapted to our notation) a well-known and central result (see e.g. Proposition 1.1 and 
Remark 1.3 in [TT]) : 

Theorem 5.3 (Definition of the WM) Consider a "semiclassical family of functions" {f^{x)}i;^(^o^i-^ 
satisfying the condition 

3Mo>0 : ||ri|L2(Mn)^Mo. (188) 

Then the family of the semiclassical WTs {M^^[/^](a;, A;)}ee(o,i) has weak-* accumulation points 
as a set of functional on an appropriate space of test-functions on phase space. When the 
accumulation point is unique (equivalently, up to the extraction of a subsequence) it will be 
called the WM associated with the semiclassical family {/^(x)}eg(o,i). 



^•^This could be seen as a "signal-processing- inspired" introduction for the semiclassical regime, see also [7]- 
The semiclassical regime, as the name shows, can also be seen as a physical regime of "large" quantum systems. 



as was the original motivation of Wigner [231 [H] • 

^^Numerical examples also offer important insights in this question. Of course, in certain problems it might 
be that some other scaling is better. We only propose this as a reasonable, general-purpose starting point. 
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Usually we will consider families with a unique accumulation point 

W'[f]{x,k)^W\x,k). (189) 

An important example is given for WKB families f^{x) = A(x)e^^^^\ where (under appropri- 
ate auxiliary assumptions, see e.g. [6J) 

W'[f]{x,k)^ \A{x)\^6{k-VS{x)). (190) 

The WM (i.e. Theorem 15.31 but also concrete examples such as that of equation (11901) 
above) is the true justification for the semiclassical scaling of the WT. 

A very successful technique in semiclassical limits has been to use the WM W^{x, k) to keep 
track of the "data of the problem" (i.e. an appropriate family of observables) in an asymptotic 
problem, i.e. for e << 1. This has been successful in many cases pTf [TT] to name but some 
landmark works. However, this approach has its own limitations, such as leading to inconsistent 
/ ill-posed problems in some cases - see [5] for a recent survey. We will also point out a couple 
of other issues here - which exist even when the WM based model can be formulated and is 
well-posed: 

• Due to the interference terms, the incorporation of e-dependent corrections to a WM 
based model is virtually impossible. This introduces a rigid scheme of the information 
that can be kept track of or not. (Indeed one might say that morally, this is at the root 
of some of the problems surveyed in [5]). The SWT offers, as we wish to show, a more 
flexible way to decide "how much detail to keep" . Some more quantitative results in this 
direction have also been presented in [21[T]. 

• The fact that we have to work with a singular object (a measure supported on a low- 
dimensional manifold) introduces many analytical as well as numerical nuances. We 
believe that regularizing to a "nice" smooth density has the potential to make many 
things easier, or even possible for the first time. A concrete, quantitative result in that 
direction is Theorem 15.81 



An important fact that we need to mention here; it is well known that in the case of critical 
smoothing, the WM is preserved. We quote the following result (in adapted notation) from 

m-. 



Theorem 5.4 (Husimi has the same weak limit as Wigner) Consider a semiclassical fam- 
ily {f} with WM 

W'[f]{x,k) -^ W\x,k) (191) 

2-k\x\^ 

in L'^(M?"')-weak. Denote g^{x) = 2^e ^ , 

H'[f]{x,k)= J fix-x')g'ik-k')W'[f]ix',k')dx'dk' = ^^,^^W'[r]ix,k). (192) 
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Then 

H'[f]{x,k) -^ W'^ix.k) (193) 

in L2(R2n).y;ea^. 

By a straightforward adaptation of the same proof as in [19j the same can be seen to hold 
in the more general case of (e-independent) crl,al as well, i.e. the SWT has the same weak 
limit as the WT, 

W'[f]{x,k) -^ W\x,k) (194) 

in L2(M2n)_^gak. 

This is important, because it shows that working with the SWT is in fact not a different 
approach than the WM. As soon as we take e ^ 0, working with the WT or with the SWT 
are indistinguishable. The difference we want to build on, is that the SWT behaves drastically 
better in several respects - e.g. numerically - than the WT in the regime < e << 1. 

A useful device in working with WTs in the semiclassical regime is the asymptotic com- 
putation of WTs of WKB functions. Indeed, computations of that kind are used in [6l [3, [H] 
to provide valuable insights - a simple one being equation (I190p . We carry out the respective 
computation for the SWT: 

Consider a WKB function of the form 

f{x) = A{x)e^'''^. (195) 

Let us suppose moreover that both A and S are analytic in a neighborhood of g G M". We 
know that the SWT of /^ will be localized near {q,p := VS{q)). The following result makes 
this precise: 

Theorem 5.5 (Asymptotic computation of the SWT of a WKB function) There is a 
neighborhood Q of the point {q,p := VS{q)) such that, V (x, k) G Q, 

\A(x)\^ -^(fc-V5(x))r-— ,^^-^— ^^(fc-VS(x)) 1 

W'[f]{x, k) = ' ^^' — P ^ i+(^D^SM)- ^ ^^^1^ 



ecri 



2\2n / i+3(£3LZ)2 5(2,))2 



^ Det 



(196) 
Let us sketch the proof of the Theorem. Let us compute M^^[/^]: 

W'[f]{x, k) = ^^ / A{x' + 5)A{x' - 5)6^'"' ^ — e ^^ ^^dx'dk'dd 

(197) 
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The stationary points of the phase are real only ii p = VS{q). When p ~ VS{q), following the 
method of [21j and since A and S are analytic we can change the path of integration in order 
to catch the complex stationary points which are given by the equations: 

- 2k' + VS{x' + 6) + VS{x' -6) = (198) 

iVS(x' + 6)- tVS(x' -6)- 2^^ = (199) 

k' — k 
-2id - 2 ^ = (200) 

Let us compute everything for k + VS(x), k — k' and x — x' small. We get 

5 = i — — (201) 



n 



and 



So 



- D^S{x') f ^) - ^ = (202) 

-k + \/S{x) + D^S{x){x' -x)-{k' -k)=0 (203) 



x' -x = -^D^S{x){k' ~k) (204) 



and 



therefore: 



^k 

2 

k' -k = -[l + ^{D^S{x)Y]-\k - VS{x)) (205) 

^k 

2 

x' -x = ^D^S{x)[l + i—D^S{x)f]-\k - VS(x)) (206) 

It is easy to check the non singularity of the Hessian of the phase. 

Finally we get the result, for k ~ VS{x), that is, \6\ ~ 0, by expanding the phase around 
the critical point. 

Now we are ready to go to our results regarding the smoothed Wigner calculus and equations. 
First of all, by obvious adaptation of the respective proof, we readily see that Theorem 14. 51 scales 
as follows: 

Theorem 5.6 (Semiclassical smoothed Wigner calculus) Let f{x), g{x) E iS(]R"), L(x,k) e 
5'(M2'^) and U = L^{x,ed^) 0. Then 

W'[L'f,g]ix,k)=C'Wy,g]{x,k) (207) 



^^See Section [5] for the scaled Weyl calculus. 
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where 



C'^w'^{x, k) = 



5,TeIR" 



where w^{x,k) = W^[f,g]{x,k) for brevity. In the special case L'^{x,k) = V{x), the respective 
expression is 

/2 
y^^y^iSx-^alS^^e^^ + '-^S, k - -S)dS (209) 

5eiR" 

All the functional analytic framework we constructed for smoothed functions should be 
scaled correctly with e; for the most part this is a very predictable exercise. The guideline is 
naturally the substitution of equation fll87p . 

The following result is an elaboration which will be particularly useful in the sequel: 

Theorem 5.7 (Semiclassical estimates for smoothed Wigner distributions) Let f^{x) E 

w'{x,k) = W'[r]{x,k). (210) 



Then the following estimate holds: 



2"||rN2 



\w'{x + iy,k + tz)\<: "" /r^ e-V-i ' ^t J . (211) 



'x'^fc 



Moreover, let mi,m2^m^,mi G (N U {0})", and denote \mi + m2 + m^ + mi^\ = m. If 
\y\, \z\ ^ ^/e, we have 



\^T^^-^^^^^w^ix + ^y,k + ^z)\ ^ F(m„ m„ m„ m,) !!(] ';:^f:i' J ^^^.^^ (212) 
Proof: First we will prove equation fl21ip . The starting point is the observation that 

W^[f^]{X, K) = T,,,^x,K [W^^x, k)] = I e-^-^'^/ (x - '-^^ g (^x + ^^ dx. (213) 

This gives us the uniform in e bound 

iw^^[r](x,K)KiiriiW.. (214) 
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Now we have 



\w'{x + iy,k + iz)\ = \ J e2'^*[(^+*2')^+('=+^^)^lu)^(X, K)dXdK\ = 



X,KeR" 

ET / 2 \r2 I „2 ^^2'^ 



!^^S£m,^(^+^) 



For equation fl212p . observe that 

^ (27r)"''| r J^mi+ni2_^m3+m4g27rJ[{x+J2;)X+{fc+i2)X]^£/j^ K)dXdK\ ^ 



^ II/'^IIl2«"')(2^)' 



2^/ l!/l%NI^ 



n I |m3+m4|d , . . |m3+m4L-i 

n 2 E rr'^) ^ r(^)^ 

d=l \ 1=0 V V fc / 



2^f |y|^ + kP 



II f^|r e ^ ^ fc 

^ F(mi, ma, mg, m4) m^/'i^^'. 



e^-^"a, ., 



2|4-n Im34-m4|+n • 



(215) 



fe||2 /Q \T?i e ^ ^ fc 

li2(lRn) 1,^/1 ; j'^N^+"^|™i+™2l+"^l'"3+™4l+" (216) 

n / \mx+m2\d , . |mi+m2|d-; 



The elementary computation which allows us to pass from the third line to the fourth, is 

+00 b^ m , \ m-l +00 „ 

/ x'^e-"" -''"rfa; = ^^ E (?) (27^) / ^'e"" rfw ^ 

r=0 '^^^ /=0 ^ '' „_ ft 

2^ 

(217) 

b^ m / \ rn—l 






The assumption |?/|, |-z| ^ \/£ implies that ^, -^ ^ 1, and therefore we can pass to the 
last line (equivalently, F{'mi,m2, 1713,1714) is independent of e). 
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The proof is complete. 

Now we are ready to see how the smoothed Wigner calculus can be approximated by dif- 
ferential operators in the semiclassical regime. This is the kind of computation necessary for 
the formulation of asymptotic SWT-based models: 

Theorem 5.8 (Semiclassical finite-order approximations to the smoothed Wigner calculus) 

Let N eN, V{x) : M" ^ M, £ G (0, 1). Assume that 

(Al) V{k) has no singular support outside {0} lj. 
(A2) 3^1 > 0, ^ Ml ^ iV + 1 such that 

\V{k)\^Ci\k\-^^' y\k\^l,k^O (218) 

(A3) For an appropriata^A (finite) constant M2 = M2{n, N) ^ min{—n — 1, —3}, 3C2 > such 
that 

\V{k)\ ^C2\k\^^^ y\k\ > 1. (219) 

Moreover, consider a "semiclassical family of wavef unctions" {f^} C iS(M"') generating the 
SWTs 

w'{x,k) = W'[f']{x,k), (220) 

for which we assume that 3Mo > such that 

||rilL2(Mn)^Mo V£>0. (221) 

According to Theorem \5.6[ 

/2 
^2niSx-^alS^y^g-^^e^^ + '-^S, k - ^S)dS. (222) 



Assuming in addition that a^ ^ 2, this expression can be approximated by differential operators, 

W^[Vf^J%x,k) = 

= E (£f E (^^D (-I)'"-' E E ^-^#^ n n d^^dl^'w^{x, k)+ 

m=0 ' ' 1=0 A e (N U {0})" B e (N U {0})" d=l d'=l 

\A\=l \B\=m-l 

+r^{x,k), 

(223) 



"'^^In fact we could handle singular support away from with no big problems; we exclude it here for simplicity. 
^^This is not the full assumption for M2. See remark 3 below, and the remarks at the end of the proof. 
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where 

||r,||i^(M2n) = 0(6"^-"), (224) 

and V{x) is the potential V{x) smoothed at scale a^, 



V2 \ f -M 



■\x~x'\'^ 



V{x) = I -^ 1 J e ^^^V{x')dx' = ^^e.V{x). (225) 

Remarks: 

1. Equation (I196p shows clearly that the L°° norm of the SWT of a WKB function is of the 
order £~". This makes the estimate of equation (12241) significant, since it can be written 

as 

||rj||i^(K2n) = 0^5^ ||#^[/^]||ioo(IK2n) j. (226) 

Moreover the same estimates as (11961) can be proved to be valid for the suitably scaled 
derivatives of Vr^[/^], making the estimate (I224p sharp. 

2. A qualitative description of the result: for appropriately (but finitely in any case) smooth 
potentials, the smoothed Wigner calculus (and therefore the "scattering term" in the 
smoothed Wigner equation (I177p ) can be approximated uniformly by a differential op- 
erator in the semiclassical regime. Of course w^{x,k) and W^lVf^, f'^]{x,k) themselves 
become unbounded pointwise as e -;> 0, but still we can approximate them strongly. 

This should be compared of course to the weak approximation of the Wigner calculus 
that is the standard device for constructing asymptotic equations for Wigner measures. 
Indeed, this result is a precise quantification of the argument that "the SWT is better 
suited to keep track of the wavefield in the semiclassical regime than the WT" . 

3. A note must be made on the selection of M2: at several instances along the proof, a 
condition of the type M2 ^ sq will appear. Some of the conditions originally appear not 
in that form, but in all cases they can be satisfied by choosing M2 small enough. The 
collection of these conditions (which depend on n, N as well) is the actual assumption 
which has to be satisfied by M2. References to all the conditions are gathered in a remark 
in the end of the proof. 

Proof: The central idea of the proof is actually very simple: we break the dS integral over a 
neighbourhood of zero and its complement, we Taylor-expand w'^{x + ^-^S, k — ^S) in equation 
(12221) around {x, k) up to order A^ and keep the remainder. Then we compute the approximation 
error of the Taylor expansion, and the contribution of the dS integral away from zero. 

First of all let us fix notations on the Taylor expansion: if g : M" ^ M is a sufficiently smooth 
function, then 

m=0 
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where the remainder can be described as follows: 39 = 9{x) G (0, 1) such that 

N+l 



Rn{x) 



n 
d=l 



(iV + 1)! 



-g[Ux) 



Now define 



g{S)=W^{x+^S,k-'-S). 



It is clear that the Taylor Theorem is applicable lj. Observe moreover that 
dsAS) 



%ea 






w%x+'-^S,k-'-S). 



Now we have the Taylor expansion 



N 



9{S) = E lS^9s. + -- + S^9s„^ ,(0) + R^^s), 



m=0 



m\ 



where, for each S, the remainder is given by 

N+l 



Rn{S) 



-N+l 



(Ar+l)!2^+i 



X; iolSdd^^ - Sddka 



w'[x + e'-^,k-ef 



(228) 
(229) 

(230) 



(231) 



^+1 / " 

,iv+i ,5-^ ^jv+1 



{N+l)\2'^^ 



1=0 



T E rt') E ^^iSdd.d E -s,dk. 



d=l 



d=l 



N+l-l 



w'{x + e'^,k-ef 



N+l 



/ 



1=0 



E ¥ n s^: 



Ae (Nujo})" 
V 1^1 = ' 



d=i 



( 



(N+i-l\ n 
\ B J Y\ qBaf)Ba 
N+l^l\ 11 ^d ^i 
Be{NU{0})" d=l 



E 
y iBi = N+i~i 



w^{x + e'^,k-ef 



N+l 



2N+1 



( 



1=0 



^daA^ 



E iU ^ra 

A e (N U {0})" d=l 

V 1^1 = ' 



/ 

n 

E in ^f C 

Be(NU{0})" d=l 

\ \B\=N + l-l 



wUx + e'^,k-e 



(232) 



-'^^Indeed, g{S) can be seen is actually an entire function; observe however that we are only interested in 
S e M", so we can use the Taylor expansion for real functions. 
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It follows therefore that, if I*?! ^ r, 

/ 



|fl„(S)K<i™^"fa? 



/=0 



Ae (Nu{o})" 






\ 



sup 

\S\^r 

A' e (Nu{0})2" 

\A'\ =N + 1 



^ B\ 

B 6 (NU{0})" 
\\B\=N + l~l J 

iecrlS 



ie...„.,....„*Ma; + ^^,^-^^ 



< 



Z=0 



sup 

\S\ < r 

A' e (N u {o})2" 
|A'| =Ar + l 



l^xi...x„fci...fc„^'^ ( ^ "^ 2^1 ^ 



e5 
2 



{ne\S\) 



N+l 



2(N+l)^ 



{N+l) 



Yy|— maxjl, (Tx I sup 

\S\^r 

A' e (N u {o})2" 

l^'l = Af + l 






\K...Xr.k,...k-.U''{X + 



£ I ^ t i^^l^ A; — — 



2 ' 



At this point we need to use Lemma IHTl To do that, we have to check that 



sal 



-f\S\ ^ V~e. 



(233) 



(234) 



One thing we will do is use the assumption cr^ ^ 2; moreover, (for reasons that will become 

1 I 1 3 I 1 

more clear below), we will set r = e '' 2{Mj+Ty^ g^j^^ therefore e\S\ ^ e" 2{a/2+i) _ gg using finally 
the constraint 

M2 ^ -3, (235) 

it follows that er ^ £2 and equation (12341) holds. In particular, the assumptions of Theorem 
15 .71 are satisfied. 

Now, using Theorem 15. 7[ and more precisely equation (12121) . as well as equation (I22ip . it 
follows that there is a constant C = C{N, n, a^, af) such that 



sup 

\S\ < r 

A' e (NU{0})2" 
\A'\ = Af+l 



l^xi...a;„fci...fc„^'' ( ^ + 2 



2 J^^ e'^+n 



(236) 



So now it follows that the remainder of the Taylor expansion of 

2 
giS) = w%x+^S,k-'-S) 



(237) 
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around S = (and for S in any case not larger than l^l ^ e 



4^2{M2 + l) 



is dominated by 



|i?;v(^)| ^//£^^"|S| 



^^^-n.| C|Af+l„^^|5|2 



e 2 



(238) 



for some constant H = H{N, n, a^, af). 

The next part of the proof is simple (if a little tedious): we break the integral of equation 
(12221) to Ji = / and I2 = J ■ We use the A^-order Taylor expansion of 



\S\^e 



"3+2(M2 + l) 



|S|>e ^+2{M2 + l) 



g{S) in Ji, and bound the error using equation (I238p . For the contribution of I2, we will use 
the estimate of equation (12 lip . Of course some more auxiliary assumptions (described in the 
statement) will come up along the way. (See also the remarks at the end of the proof). 
The first contribution to the error comes from 



Ei= / e'"*^"-^'^'^'y(5)i?,v(5)rf5. 

Using the previous results and the assumptions of the Theorem we have 

\Ei\^He^-^ f \V{S)\\S\''+^dS ^ 



<Hs'^- 






-Mi+N+l 



dS + C2 



I 



\S\M2+N+l^g 



He 2 



r(tT 



r=0 



11 1 
3^2(M2 + l) 



-Mi+N+n 



dS + C2 J r^''+''+''dS 



r=\ 



Here we make use of the assumption Mi ^ A^ + n, and moreover we assume 

M2 + A^ + n + 1 ^ 0. 
Now we are able to proceed to 



\E,\^He^- 2^" 



r(f) 



Ci 



A^ + ra + 1 - Ml 



+ C2 



( tiU I'^MtiM'i/' 1| 1 \ 




^\l\i2 1 iV in 1 ij(^ ^1 2(M2 + l)j 


- 1 


M2 + A^ + n + 1 





Clearly the best we can ask for here is 

(M2 + AT + n + 1) 



'4 2(M2 + 1) 



>0, 



(239) 



(240) 



(241) 



(242) 



(243) 
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N+1 



SO that the total bound for Ei is controlled by e^ " (it is obvious that for M2 small enough 
the inequality fl243p holds). 

On the contribution of I2, making use of Theorem 15.71 we observe that 



l/.l ^ 



/ \V{S)\e-^^^^\w{x+'^S,k-'^S)\dS ^ 



\S\>e * + ^p4+T7 



(244) 



< 



C2Mg 27rt 



s r 



/ 



r 



M2+n-l 



dr 



2/ 1, 1 

r=s 3+2(Af2+l) 



We will have to assume 



Mo ^ -l-n 



(245) 

for the integral to exist (observe however that this is not automatically enough for lim J2 = 0); 
with that we get 



\TJ < '^2Mg 27r7 2 4 l-2(M2 + l) 

M2I ^ (^M2-n)r(f)^ • 



(246) 



JV+1 



Like earlier, we ask not only that lim/2 = 0, but that I/2I is controlled by e 2 ". This amounts 

to 

_l_M^ n-1 jV+l_„. ,24,) 

2 4 2(M2 + 1) 2 ^ ^ 

Now observe that equation (12321) gives the structure of all the terms in the Taylor expansion, 
not only the remainder. That is, the order-m term of the Taylor expansion (I23ip is given by 



2\^ I 1 \m-l 



TUS) = 9E i^^l) (-1 



1=0 



/ \ 




/ \ 


yle(Nu{0})" <i=i 
V 1^1 = ' / 




2^ B! 11 -^d % 

B6(NU{0})" d=l 
\ |_B| = m - / / 



w^ (x, A;) . 



(248) 



So what we have shown so far is that 






^2TTiSx-^C7iS' 



N 



vis)z 



m=0 



m!2" 



d=l 



dSw'^{x, k) + r^{x, k) 



N rri 



EPE(^^^)'(-i 



2V ( 1 \m-l 



E 



E ^^^ n n di^d^/w^ix, k) + f,(x, A;) 



m=0 i=0 



A6(NU{0})" B6(NU{0})'^ 
|A|=« \B\=m-l 



d=l d'=l 



(249) 
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where the coefficients D{A) are given by 

/n 
e-^'''^''-'T-'^'v{S)Y[S^dS, (250) 

J 1 



|5Ke 3-^2{A.f2+l) 

and 



1^ 1 d=l 



JV+l 



\\fs\\L^iR^«) = 0(^6 — -"j. (251) 

Observe that the coefficients D{A + B) are truncated versions of the derivatives of the 
smoothed potential. Indeed, denote 



_ n 

D{A){x) = J e"2-^5"-'?"'^'y(5) n Sf'dS 

d=l 



1,2 



(252) 



then 






h{A){x) = D{A){X) - D{A){X) = f e-^niSx-^aiS^y^g^ ]-J gA.^g^ 



_i+ 1 d=l 

\S\>6 5+2(M2+l) 



(253) 



This difference is small, and the final part of the proof consists in showing that we can substitute 
D{A) for D{A) and have a similar error estimate as in equations (12491) . (I25ip . To that end, it 
suffices to show that 3^ > such that Vx e M" and V^ G (N U {0})", |^| ^ A^, 

\h{A){x)\^Ce'^. (254) 

Then (using Theorem 15.71 once more) the error introduced by substituting D{A) for D{A) in 
equation (12491) will be dominated by 

N m I r ^ ^ nM "■ "■ D 

E f^EK') E E ^^^^ n n \di^d^,fw^{x,k)\ + \ux,k)\ ^ 

m=0 1=0 A6(NU{0})" B6(NU{0})" d=ld'=l 

\A\ = l \B\=m-l 

N , 

m=0 ^ 

(255) 
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So let us prove equation fl254p : 



l3{A)ix) ^ 
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e-'f-i 


^"\v{s) n 


^.1^' 


*(i5^ 




\S\>e' 


1 1 1 

3 + 2{M2 + l) 




d=l 








<C2 
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^-^.lS^^gM2+n\A\^g^ 






\S\>6~ 


1 1 1 
3"^2(Af2 + l) 










= 


r=e 


+00 

/ 

1 1 1 


g 2 ^ 


r2^M2+n|A|+n 


-^dr 
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n 

r(f)(- 


C2 

-e7r<72)^'^2+n|A|4 


-n-1 

1 


+00 

/ 

+ 2(M2 + 1) 


oM2+n|A|+n- 


1 _£2! 

e 2 


alr^^^ ._ 
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27rT 

"r(t) 


C2 p.M2+n\A\+ri- 


2 _£2:„2 2 

g 2 a: 


1 I 

"1+2C 


1 ^ 

M2+1) ^ 



(256) 



< C" ^-M2-n\A\-n+l 



TT I 7ro-?e^'*"2(M2 + i) 

-nM2+n|A|+n-2 



e 2 



eir 2_7^ Ai9 + 1 



■<t:££^ *^2 



where of course Hs{x) is the Hermite polynomial of order s. By strengthening the assumption 
M2 ^ —3 (which already has appeared in equation (I235p ) to 



M2 < -3 



(257) 



o I 1 

we get that e* 2(M2+i) = o(l), and therefore only the zero order term of H M2+n\A\+n-2 has to 
be considered in the last line of equation (I256p . Observe moreover that 



s^aje^ ^^^i+T <- ]_ 



g 2 ^ 



Using these observations, equation (12561) implies 

\h{A){x)\ ^ C^£-A^2-n|A|-„+l ^ (jn^-M2-nN-n+l (258) 

(remember that \A\ ^ A^). Now asking that equation (12541) holds is equivalent to asking 



e'M2-nN-n+l ^e^ ^ _M2 - nN - u + I ^ ^ ^ 



^ M2^-N{n+ i) 



(259) 



n 



2' 



The proof is complete. 
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Remarks on the choice of the parameters: 

The constraints for M2 come up during the proof in equations (12351) . (12411) . (I243p . (I247p . 
(1257p . (I259p . Clearly, each of these constraints can be satisfied for M2 small enough, depending 
on A^, n. 

The constraint for a^ came up right after equation (I234p . 

The constraint for Mi appeared right after equation (I240p . 

Theorem 15.81 result must be compared with its counterpart for the WT, quoted here (in 
adapted notation) from [TT] : 

Theorem 5.9 (Semiclassical finite-order approximations to the Wigner calculus) Let 

p{x, k) e C°°(M2") satisfy, for some M( > 

V aeN2"|a,"^_„,^,.,,„p(x,fc)| < Ca{l + \k\Y'K (260) 

Assume also ||/^||l2(R"), ||fi'^||L2(iR") ^ M) Ve > 0. Then 

W'[p{x,ed,)f',g%x,k) = p{x,k)W'[P,g%x,k) + 

(261) 
+ iri E [dkA^. k)d,^W^[f^, g%x, k) - d^^x, h)d,^Wy^, g%x, k)] + eV,, 

d=l 

where r^ is hounded in iS'(M^") as e ^ 0. 
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